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ABSTRACT 


Certain  ORION  missions  may  require  three  axis  stabilization.  Since  ORION'S 
physical  size  severely  limits  its  onboard  fuel  storage  capability,  passive  stabilization 
techniques  warrant  investigation.  This  paper  shows  the  development  of  linearized 
equations  of  motion  and  regions  of  stability  with  respect  to  gravity  gradient  stabiliza¬ 
tion.  Gravity  gradient  stabilization  by  itself  provides  little  yaw  restoring  torque;  there¬ 
fore,  additional  torque  generating  devices  are  necessary  to  augment  the  gravity  gradient 
effect.  Control  moment  gyros,  reaction  wheels,  and  magnetic  torquers  will  be  investi¬ 
gated  as  to  their  suitability  for  ORION. 


Accession  Fo 

NTIS  CFA*I 
DTIC  XAfi 
Unannounced 

J.i3l  if  tout  to 

- - — • 

P  / 

ur 

□ 

□ 

|  o  > 

Distribution/ 

1  Availability  Codes 

t 

101st 

r' 

Avail 

Spec 

and/or 

lal 

TABLE  OF  CONTENTS 


I.  INTRODUCTION  . 1 

A.  BACKGROUND  ON  ORION  . 1 

1.  Background  on  attitude  control  . 2 

II.  GRAVITY  GRADIENT  STABILIZATION  THEORY  . 5 

A.  COORDINATE  DEFINITION  . 5 

B.  COORDINATE  TRANSFORMATION  . 6 

1.  Direction  cosines  method  . 6 

2.  Quaternions  . 7 

3.  Euler  angles  method  . 9 

C.  RIGID  BODY  EQUATION'S  . 13 

1.  Derivative  of  a  Vector  in  a  Rotating  Coordinate  System  . 13 

2.  General  Particle  Motion  . 14 

3.  Momentum  of  a  Rigid  Body  . 15 

4.  Equations  of  motion;  Euler's  Moment  Equations  . 17 

D.  GRAVITY  GRADIENT  EFFECT  . 20 

E.  DESIGN  FOR  STABILITY  . 23 

1.  Stability  and  Frequency  Analysis . 23 

III.  ENVIRONMENTAL  DISTURBANCES  AND  RESTORING  TORQUES  .  .  27 

A.  ENVIRONMENTAL  DISTURBANCES/TORQUES  . 27 

1.  Aerodynamic  disturbance  . 27 

a.  Drag  . 27 

b.  Aerodynamic  torque  . 27 

2.  Rigid  Body  Solar  Pressure  Torque  . 30 

B.  RESTORING  TORQUES  . 31 

1.  Thrusters  . 31 

2.  Momentum  Exchange  Devices  . 32 

a.  Reaction  wheels  . 33 

3.  CMG(s)  . 35 

4.  Magnetic  torquers  . 35 


1 


IV.  CHAPTER  4  . 39 

A.  GRAVITY  GRADIENT  STABILIZATION  RESULTS  . 39 

1.  Attitude  Control  . 39 

2.  Various  results  based  on  different  configurations  . 39 

a.  Effects  of  extending  booms  . 39 

b.  Gravity  gradient  effect  alone  . 39 

c.  Effects  of  Reaction  VVheel(s)  . 42 

d.  Effects  of  Magnetic  torquers  . 5S 

3.  Gravity  Gradient  capture  {from  tumbling  mode)  . 79 

V.  CHAPTER  5  . 81 

A.  CONCLUSIONS  . 81 

1.  Gravity  gradient  stabilization  (only)  . SI 

2.  Gravity  gradient  with  3  reaction  wheels  . SI 

a.  Advantages  . 81 

b.  Disadvantages  . 81 

3.  Control  moment  gyro  (CMG)  . 81 

a.  Advantages  . 81 

b.  Disadvantages  . 81 

4.  Gravity  gradient  with  three  magnetic  torquers  . 82 

a.  Advantages  . 82 

b.  Disadvantages  . 82 

5.  Thrusters  . 82 

a.  Advantages  . 82 

b.  Disadvantages  . 82 

6.  Overview  . 82 

APPENDIX  GRADIENT  STABILIZATION  PROGRAM  . 84 

LIST  OF  REFERENCES  . 89 

INITIAL  DISTRIBUTION  LIST  . 90 


v 


LIST  OF  FIGURES 


Figure  1.  Nasa  Get-Away-Special  (GAS)  Cannister  [Ref.  1,  p.  30]  . 1 

Figure  2.  Cross  section  of  GAS  Cannister  [Ref.  1,  p.  42]  . 2 

Figure  3.  Summary  of  typical  Xavy/STP  payload  requirements  . 3 

Figure  4.  Passive  Stabilization  Methods  [Ref.  4,  p.  19]  4 

Figure  5.  Active  Stabilization  Methods  [Ref.  4,  p.  19]  . 4 

Figure  6.  Coordinate  Definition  for  Gravity  Gradient  Systems  . 5 

Figure  7.  Coordinate  Transformation  Methods  [Ref.  4,  p.  412]  . 6 

Figure  8.  Rotation  about  the  Euler  Axis  . 7 

Figure  9.  Euler  Angles  . 10 

Figure  10.  Rotation  about  the  Z  axis  through  the  Angle  . 10 

Figure  11.  Table  of  Euler  angle  combinations  [Ref.  4,  p.  764]  . 12 

Figure  12.  Vector  in  a  rotating  system  [Ref.  6,  p.  107]  13 

Figure  13.  Translational  motion  of  a  particle  [Ref.  6,  p.  107] . 14 

Figure  14.  Gravity  gradient  torque  [Ref.  6,  p.  113]  . 21 

Figure  15.  Gravity  Gradient  torque  Magnitudes  vs  Altitude  [Ref.  3,  p.  5S]  . 23 

Figure  16.  Gravity  Gradient  Stability  Regions  . 26 

Figure  17.  Environmental  Disturbance  Torques  [Ref.  4,  p.  17]  27 

Figure  18.  Aerodynamic  effects  for  cylindrical  satellite  [Ref.  3,  p.  5S]  . 2S 

Figure  19.  Aerodynamic  pressure  as  a  function  of  altitude  [Ref.  8,  p.  456]  . 29 

Figure  20.  Aerodynamic  torque  for  a  1  inch  offset  as  a  function  of  altitude . 30 

Figure  21.  Rigid  body  solar  pressure  torques  [Ref.  3,  p.  61] . 31 

Figure  22.  Fuel  required  for  orbit  transfer  based  on  a  2501b  satellite.  [Ref  2,  p.  23]  32 

Figure  23.  3  axis  reaction  wheel  control  scheme  [Ref.  6,  p.  149]  34 

Figure  24.  Twin  gyro  controller  with  one  drive  motor  [Ref.  8,  p.  415]  . 35 

Figure  25.  Magnetic  Moment  Due  to  a  Current  Loop  [Ref.  4,  p.  204]  . 36 

Figure  26.  Some  existing  magnetic  moment  schemes  [Ref.  3,  p.  205]  . 37 

Figure  27.  Earth's  magnetic  field  modeled  as  a  tilted  dipole  [Ref.  6,  p.  148]  . 37 

Figure  28.  Effect  of  extending  boom  on  moments  of  inertia  . 40 

Figure  29.  Gravity  gradient  effect  only  Ix=  15,  Iy=  18,  Iz=  3  . 41 

Figure  30.  Gravity  gradient  effect  only  Ix=  95,  Iy=  98,  Iz=  3  . 42 

Figure  31.  Teldrix  Reaction  Wheel  Technical  Data  . 43 


VI 


Figure  32.  Bendix  Aerospace  Reaction  Wheel  Technical  Data . 44 

Figure  33.  Gravity  gradient  efTect  with  3  RW  (X  axis)  lx  =  15.  Iy=  18.  Iz=  3  ....  45 

Figure  34.  Gravity  gradient  effect  with  3  RW  (V  axis)  Ix=  15.  Iy=  18.  Iz=  3 . 46 

Figure  35.  Gravity  gradient  efTect  with  3  RW  (Z  axis)  lx  =  15,  ly  =  18.  Iz=  3  . 47 

Figure  36.  Gravity  gradient  efTect  with  3  RW  (X  axis)  lx  =  95,  I y  =  9 S ,  Iz=  3  ....  48 

Figure  37.  Gravity  gradient  efTect  with  3  RW  (V  axis)  lx  =  95,  Iy=  98,  Iz=  3 . 49 

Figure  38.  Gravity  gradient  efTect  with  3  RW  (Z  axis)  lx  =  95,  Iy=  98,  Iz=  3  . 50 

Figure  39.  Gravity  gradient  efTect  with  3  RW  (X  axis)  Ix=  15,  Iy  =  IS,  Iz=  3  ....  52 

Figure  40.  Gravity  gradient  efTect  with  3  RW  (Y  axis)  lx=  15,  ly=  18,  lz=  3 . 53 

Figure  41.  Gravity  gradient  efTect  with  3  RW  (Z  axis)  Ix=  15.  Iy=  IS,  lz=  3  . 54 

Figure  42.  Gravity  gradient  efTect  with  3  RW  (X  axis)  lx  =  95,  Iy  =  98,  Iz  =  3  ....  55 

Figure  43.  Gravity  gradient  efTect  with  3  (Y  axis)  RW  Ix  =  95,  I y  =  98,  Iz  =  3  . 56 

Figure  44.  Gravity  gradient  efTect  with  3  RW  (Z  axis)  Ix=95,  Iy=9S,  Iz=  3  . 57 

Figure  45.  Gravity  gradient  with  magnetic  torquing  (X  axis)  . 60 

Figure  46.  Gravity  gradient  with  magnetic  torquing  (Y  axis)  . 61 

Figure  47.  Gravity  gradient  with  magnetic  torquing  (Z  axis)  . 62 

Figure  48.  Gravity  gradient  with  magnetic  torquing  (X  axis)  . 63 

Figure  49.  Gravity  gradient  with  magnetic  torquing  (Y  axis)  . 64 

Figure  50.  Gravity  gradient  with  magnetic  torquing  (Z  axis)  . 65 

Figure  51.  Gravity  gradient  with  magnetic  torquing  (X  axis)  . 67 

Figure  52.  Gravity  gradient  with  magnetic  torquing  (Y  axis)  . 68 

Figure  53.  Gravity  gradient  with  magnetic  torquing  (Z  axis)  . 69 

Figure  54.  Gravity  gradient  with  magnetic  torquing  (X  axis)  . 70 

Figure  55.  Gravity  gradient  with  magnetic  torquing  (Y  axis)  . 71 

Figure  56.  Gravity  gradient  with  magnetic  torquing  (Z  axis)  . 72 

Figure  57.  Gravity  gradient  with  magnetic  torquing  (X  axis)  . 73 

Figure  58.  Gravity  gradient  with  magnetic  torquing  (Y  axis)  . 74 

Figure  59.  Gravity  gradient  with  magnetic  torquing  (Z  axis)  . 75 

Figure  60.  Gravity  gradient  with  magnetic  torquing  (X  AXIS)  . 76 

Figure  61.  Gravity  gradient  with  magnetic  torquing  (Y  AXIS)  . 77 

Figure  62.  Gravity  gradient  with  magnetic  torquing  (Z  AXIS)  . 78 

Figure  63.  Gravity  gradient  capture  by  extension  of  boom  . 80 


vii 


I.  INTRODUCTION 


A.  BACKGROUND  ON  ORION 

The  philosophy  behind  the  design  and  development  of  the  Orion  satellite  is  being 
driven  in  part  by  economic  considerations.  Small,  affordable,  general  purpose  satellites 
must  be  developed  without  sacrificing  reliability.  The  preliminary  groundwork  on  the 
feasibility  of  the  ORION  concept  and  design  is  detailed  in  [Ref.  1]. 

Orion  is  currently  being  designed  to  be  launched  from  a  Get  Away  Special  (GAS) 
cannister  aboard  the  Space  Shuttle.  The  GAS  cannister  physically  limits  the  size  of 
ORION  to  approximately  35"  in  height  and  19"  in  diameter.  Figure  1  shows  the  Get- 
Away  Special  Cannister  attached  to  the  side  of  the  shuttle  bay. 


Figure  1.  NASA  Get-Away-Special  (GAS)  Cannister  [Ref.  1,  p.  301 


The  Get-Away-Special  (GAS)  Cannister  is  designed  to  launch  its  satellite 
autonomously  on  command.  A  cross  sectional  view  of  the  (GAS)  Cannister  is  shown  in 
Figure  2  on  page  2. 

ORION  should  be  able  to  meet  the  requirements  of  standard  or  typical  payloads. 
Some  of  these  requirements  are  listed  by  Reference  2  and  given  in  Figure  3  on  page  3. 
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Figure  2.  Cross  section  of  GAS  Cannister  [Ref.  1,  p.  42j 
1.  Background  on  attitude  control 

The  success  of  any  satellite  is  dependent  on  its  attitude  control  system.  To  be 
mission  capable,  the  vehicle  must  first  acquire  a  desired  or  specified  orientation  in  space, 
and  then  maintain  this  orientaion  within  mission  limits  given  such  external  factors  as 
magnetic  fields,  aerodynamic  drag,  solar  pressure,  gravitational  effects,  and  other  dis¬ 
turbances.  Deviations  from  desired  orientations  are  detected  through  the  use  of  sensors. 

Considerations  for  an  attitude  control  system  include: 


(1) 

pointing  accuracy 

(2) 

fuel  consumption 

(3) 

thermal 

(4) 

power  constraints 

- 

(5) 

satellite  lifetime 

(6) 

response  time 

Mass  32  lbm 

Volume  2.36  ft3 

Power  34  watts 

DataRate  5000  bits/sec 

Orbit  200-800  nm  circular 

Inclination  o°-300  or  60*-120* 

Instruments  Particle  counter  or  Lens 

Figure  3.  Summary  of  typical  Navy/STP  payload  requirements 

I  he  weight  given  to  each  of  these  considerations  ultimately  lies  in  the  specific  mission 
requirements  designated  for  each  satellite.  The  direction  of  this  thesis  will  be  to  inves¬ 
tigate  the  feasibility  of  achieving  sufficient  performance  utilizing  a  gravity  gradient  con¬ 
trol  system,  lhc  standard  of  +/-  1  degree  for  each  axis  will  be  used.  This  criteria  is 
sufficient  to  enable  ORION  to  satisfy  most  STP  (space  test  program)  mission 
requirements  ! Ref.  1,  p.  93J 

To  achieve  stabilization,  many  methods  of  generating  controlling  torques  have 
been  developed.  These  methods  are  divided  into  either  active  or  passive  categories.  Each 
of  these  categories  has  their  relative  advantages  and  disadvantages  [Figure  4  on  page  4 
and  Figure  5  on  page  4). 

In  this  study,  the  primary  focus  is  on  a  gravity  gradient  stabilization  augmented 
by  other  torque  generating  devices  to  achieve  the  desired  +  /- 1  degree  goal  for  each  axis. 
A  discussion  of  the  equations  of  motion  and  gravity  gradient  stabilization  theory  is 
presented  in  Chapter  2.  Disturbance  torques  and  other  torque  generating  devices  are 
discussed  in  Chapter  3.  Chapter  4  presents  the  results  of  stabilization  schemes  derived 
in  Chapters  2  and  3.  Observations  and  conclusions  are  discussed  in  Chapter  5. 
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SOLAR  RAOIAIION  STABILIZED 


AtnOOVNAMiC  STABILIZED 


MAONf  TIC  STABILIZED  WITH 
PERMANENT  MAGNET 


ADVANTAGES 

DISADVANTAGES 

Simplf  Frrfcnvf  n f a n i  y  anywhere  in  any 
ORIENTATION.  MAINTAINS  ORif  NIAtlON  IN  INERTIAL 

space 

CENTRIFUGAL  fORCE  REOUIRES  RTnuCTlIRAL 
STABILITY  ANO  SOME  RlC.IOITY  SENSORS 

ANO  ANTENNAS  CANNOT  GENERALLY  REMAIN 
POINTED  AT  A  SPECIFIC  INERTIAL  TARGET 
WOBBLE  INUTATIONI  IF  NQT  PROPERLY  BAL 
ANCED  ORift  DUE  TO  ENVIRONMENTAL 
TORQUES 

MAINTAINS  STABLE  ORIENTATION  RELATIVE  TO 
CENTRAL  BOOY  NOT  SUBJECT  TQ  DFCA  Y  OR  DR  IT  T 

DUE  TO  ENVIRONMENTAL  TORQUES  UNLESS  ENVIRON 
MENT  CHANGES 

LIMITED  TO  1  OR  7  POSSIBLE  OR»EN*  ATionS  EF 
FECTIVf  ONLY  NEAR  MASStVF  CENTRAL  BOOY 
IE  G  EARTH  MOON.  ETC  1.  REQUIRES  LONG 
BOOMS  OR  CLONGATEO  MASS  O'STR'BUTION 
SUBJECT  TO  WOBBLE  ILIBRATioni  CONTROL 
LIMITFO  TO  -l  OEG  -PROBLEM  OF  THERMAL 
GRADIENTS  ACROSS  BOOM 

CONVENIENT  TOR  POWFR  GENERATION  By  SOLAR 

CELLS  OR  SOLAR  STUOIES 

LIMITED  TO  HiGH-AL  TIT  UO  E  OR  IN  T  EnPL  anE 
TARY  ORBITS.  LIMITED  ORIENTATIONS  AL 
LOWED 

J  SPECIAL  PURPOSE  METHODS  -  HIGHLY  MISSION  AND  STRUCTURE  DEPENDENT  l 
\  IN  ALL  OF  IMfIR  CHARACTERISTICS  | 

Figure  4.  Passive  Stabilization  Methods  [Ref.  4,  p.  19] 


MAGNETIC  (ELECTROMAGNETS! 


REACTION  WHEELS* 


ALTERNATIVE  THRUSTERS  »ON 
OR  ELECTRIC 


ACTIVE  SOLAR.  AERODYNAMIC. 
OR  gravity  gradient 


FLEXIBLE  AND  EAST.  USED  IN  ANY  ENVIRONMENT; 

POWERFUL 


USUALLY  LOW  POWER  MAY  BE  DONE  WITHOUT  USING 
CONSUMABLES  BY  USE  OF  SOL  An  POWER 


PARTICULARLY  OOOO  FOR  VARIABLE  SPIN  RATE  CON 
TROL  FAST.  FLEXIBLE.  PRECISE  ATTITUDE  CONTROL 
ano/qr  stabilization 


USES  CONSUMABLE  IFUILI  WITH  LIMITED 
SUPPLY  AVAILABLE  DUE  TO  FUEL  WEIGHT.  TOO 
POWERFUL  FOR  SOME  APPLICATIONS  (. «  . 
RELATIVELY  COARSE  CON  TROL  I  COMPLEX  AND 
EXPENSIVE  PLUMBING  SUBJECT  TO  FAILURE 

SLON  NEAR  EARTMONLY  APPLICABILITY 
UMITEO  BY  DIRECTION  OF  THE  EXTERNAL 
MAGNETIC  FIELO:  COARSE  CONTROL  ONL  < 
(BECAUSE  OF  MAGNETIC  FIELD  MOOCL  UN 
CERTAINTIES  AND  LONG  TIME  CONSTANTS! 


REQUIRES  RAPIDLY  MOVING  PARTS  WHICH 
IMPLIES  PROBLEMS  OF  SUPPORT  ANO  FRIC 
TION.  MAY  NEED  SECOND  CONTROL  SYSTEM 
TO  CONTROL  OVERALL  ANGULAR  MOMENTUM 
I  MOMENTUM  OUMPING  I  IN  RESPONSE  TO 
CUMULATIVE  CHANGES  BY  ENVIRONMENTAL 
TORQUES  EXPENSIVE 


PRIMARILY  JPFC'Al  PURPOSE -LESS  CxrfR'.ENCE  WITH  THESE  THAN  WITH  THOSE  | 
LISTED  ABOVE  CHARACTERISTICS  ARE  HIGHLY  MISSION  OEPENQENT;  SOME 
SYSTEMS  MAY  SEE  MORE  USE  IN  THE  FUTURE  AS  FURTHER  EXPERIENCE  IS  GAINED 


•REFERS  TO  ANY  DEVICE  THAT  MAY  BE  USED  IN  A  PROCESS  TO  EXCHANGE  ANGULAR  MOMENTUM  WITH  THE  SPACECRAFT  BOOY 


Figure  5.  Active  Stabilization  Methods  [Ref.  4,  p.  19] 


II.  GRAVITY  GRADIENT  STABILIZATION  THEORY 


A.  COORDINATE  DEFINITION 

Before  beginning  any  discussion  on  attitude  control  a  means  of  defining  the  coordi¬ 
nate  system  is  necessary.  The  standardized  coordinate  axes  for  gravity  gradient  systems 
arc  as  follows:  the  x  axis  points  in  the  direction  of  satellite  motion,  the  v  axis  is  normal 
to  the  orbit  plane,  and  the  z  axis  points  towards  earth's  geocenter.  These  axes  corre¬ 
spond  to  the  roll,  pitch  and  yaw  axes,  respectively.  Figure  6  from  Reference  3  defines 
roll  (</>),  pitch  (6).  and  yaw  (i [/)  attitude  errors  as  well  as  showing  the  standardized  co¬ 
ordinate  system. 
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B.  COORDINATE  TRANSFORMATION 

The  relationship  between  two  coordinate  systems  is  essential  for  attitude  control. 
There  are  methods  available  to  describe  this  relationship  of  which  directional  cosines, 
quaternions,  and  Euler  angles  will  be  briefly  discussed.  Figure  7  on  page  6  shows  various 
methods  of  coordinate  transformations  along  with  their  respective  advantages  and  dis¬ 
advantages. 


advantages 

DISADVANTAGES 

COMMON  APPLICATIONS 

direction 

COS'Nt 

MATRIX 

A  -  IA..I 

NO  SINGULARITIES 

NO  TRIGONOMI  THIC  FUNCTIONS 
CONVENIENT  PROOUCT  RULE 

F  OR  SUCCESSIVE  ROTATIONS 

SIX  REOUNOANT  PARAMETERS 

IN  ANALYSIS  TO  TRANSFORM 
VECTORS  FROM  ONE  REFER. 

INCE  FRAME  TO  ANOTHER 

EULER 

AXI5/ANGLC 

c  ♦ 

Clear  physical  interpretation 

ONE  REOUNOANT  PARAMETER 
AXIS  UNOEf  INEO  WHEN  SIN  *-© 
TRIGONOMETRIC  FUNCTIONS 

COMMANDING  SLEW  MANEUVERS 

IUIIR 

SYMMETRIC 

P  A  RAM  E  TEAS 
IQUATERNIONI 

M 

no  SINGULARITIES 

NO  TRIGONOMETRIC  FUNCTIONS 
CONVENIENT  PRODUCT  RULE  FOR 
SUCCESSIVE  ROTATIONS 

ONE  RCOUNDANT  PARAMETER 

NO  OBVIOUS  PHYSICAL  INTER¬ 
PRETATION 

ONSOARO  inertial  navigation 

ones 

VECTOR 

• 

NO  REOUNOANT  PARAMETERS 

NO  TRIGONOMETRIC  FUNCTIONS 
CONVENIENT  PRODUCT  RULE  FOR 
SUCCESSIVE  ROTATIONS 

INFINITE  FOR  ItO-OCG  ROTA¬ 
TION 

ANALYTIC  STUDIES 

EULER 

ANGLES 

#.  f,  * 

NO  REOUNOANT  PARAMETERS 
PHYSICAL  INTERPRETATION  IS 

CLEAR  IN  SOME  CASES 

TRIGONOMETRIC  FUNCTIONS 
SINGULARITY  AT  SOME  $ 

NO  CONVENIENT  PROOUCT  RULE 
FOR  SUCCESSIVE  ROTATIONS 

ANALYTIC  STUOICS 

INPUT /OUT  PUT 

ONBOARD  ATTtTUOf  CONTROL  OF 

3-AXtS  stabilized  spacecraft 

Figure  7.  Coordinate  Transformation  Methods  [Ref.  4,  p.  412] 


1.  Direction  cosines  method 

7  he  directional  cosine  method  involves  a  3  x  3  rotation  matrix  A()  where 


^12 

•^13 

X'x 

X>y 

X'i 

AV- 

^2. 

A22 

^23 

= 

r.x 

Y-y 

Y  •  2 

^32 

1 

ro 

rn 

Z  •  x 

Z-y 

Z  •  z 

where 


X  •  .r  =  cos  9  cos  ip 

X  »y  =  cos  ip  sin  9  sin  0  +  sin  ip  cos  <p 
X  •  i  =  sin  \p  sin  <p  -  cos  ip  sin  9  cos  0 

Y  *x  —  —  sin  ip  cos  9 

Y  »y  =  cos  4>  cos  ip  —  sin  ip  sin  9  sin  <p 

Y  •  z  =  sin  ip  sin  9  cos  tp  4-  cos  ip  sin  0 
Z  •  x  =  sin  9 

Z  »y=  —  cos  9  sin  0 
Z  •  z  =  cos  0  cos  9 


6 


2.  Quaternions 

The  quaternion  method  or  as  it  is  sometimes  called,  the  Euler  Symmetric  Pa¬ 
rameter  method,  uses  the  Euler  axis  and  rotation  angle  (a  single  rotation  angle  a  about 

A 

a  single  axis  £)  method  to  describe  the  relationship  between  two  coordinate  systems 
(Figure  8  on  page  7). 


Figure  8.  Rotation  about  the  Euler  Axis 


The  quaternion  is  a  four  component  vector  containing  the  same  information  as  the  Euler 
angle/axis  transformation  in  which  one  coordinate  axis  is  related  to  another  and  is  of  the 
form 

qj,  q3k,  ft]  (3) 

where 

q2  =  Ey  sin  y 
siny 
?4-cosy 

The  symbols  i,  j,  and  k  satisfy  the  following  conditions  [Ref.  4,  p.  758J 


•2  .2  ,2  , 
i  =j  =  k  =1 

ij  =  -ji  =  k 

jk  =  —kj  =  i 

ki  =  -ik  =j 


and 


2  .  2  .  2  ,  2  , 

<7l  +  <?2  +  <?3  +  <74  =  1 

When  used  for  control  purposes,  [Ref  5]  defines  the  first  three  elements  of  the 
quaternion  (qu  q2,  qh)  as  the  respective  roll,  pitch,  and  yaw  rotational  errors: 

pitch >  = 

=  2^£x  sin-|-  ,  Ey  sin-^- ,  Ez  sin-y  j  (7) 

s:[a£x,  aEy,  o.E2~\ 

The  primary  advantage  of  the  quaternion  method  is  that  the  computational  time 
can  be  reduced  by  more  that  40%  over  the  equivalent  directional  cosine  matrix  solution 
in  operations  that  require  successive  coordinate  system  rotations  [Ref  5].  The  transfor¬ 
mation  from  coordinate  system  A  to  coordinate  system  B  can  be  described  by  the 
quaternion  q  =  [<7,/,  qj,  qj,  q 4],  and  the  second  transformation  from  system  B  to  C  can 
be  described  by  q'  =  [q\i,  q'j,  q'j,  The  transformation  from  system  A  to  system 
C  may  be  described  by  the  quaternion  q"  =  [ \q'\i ,  q"J,  q" 3,  q'\~\  where 

q"  =  qq'  =  (qy  i  4-  qj  +  q3k  +  qj){q' ,  1  +  q'j  +  q’3k  +  q\)  (8) 

Multiplying  out  and  substituting  in  the  conditions  stated  above,  the  matrix  form  be¬ 
comes 


<?”. 

— 1 

-c. 

1 

_ 1 

<7i 

<r  2 

~<7  3  9  4  9  1  ?  2 

<72 

<7  3 

<7  2  — <7 1  <7  4  ?  3 

<73 

q'\ 

-?'l  -<7'2  -?'3  ?'4_ 

.<74. 

Because  of  the  quaternion  property  of  interchangeability,  a  total  of  16  multiplications 
and  12  additions  are  required  to  completely  define  the  quaternion  describing  the  trans¬ 
formation  from  coordinate  system  A  to  coordinate  system  C.  To  accomplish  the  same 
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transformation  using  directional  cosines  would  require  multiplying  two  3x3  matrices, 
involving  27  multiplications  and  IS  additions.  The  computational  time  savings  for  suc¬ 
cessive  rotations  transformations  becomes  apparent.  The  major  drawback  for  the  use 
of  quaternions  is  that  the  numbers  of  the  quaternion  do  not  physically  represent  the 
transformation  from  one  system  to  another. 

3.  Euler  angles  method 

The  Euler  angles  utilize  three  different  rotations  (angles),  defining  the  orien¬ 
tation  of  a  body  with  respect  to  an  inertial  reference  frame.  These  angles  are  defined  as 
4>  (about  the  x  axis),  6  (about  the  y  axis),  and  \p  (about  the  z  axis).  These  axes,  by  es¬ 
tablished  convention,  may  also  be  referred  to  as  1  (for  x),  2  (for  y),  and  3  (for  z)  axes, 
respectively. 

The  321  Euler  transformation  described  below  starts  with  the  inertial  coordinate 
system  XYZ.  Figure  9  on  page  10. a.  shows  the  first  rotation  about  the  Z  axis  through 
an  angle  >p,  which  produces  xLv,Z.  Figure  9  on  page  lO.b.  shows  the  rotation  of  6  about 
y,  ,  producing  xy,r,.  Figure  9  on  page  lO.c.,  showing  the  rotation  of  </>  about  x.  produces 
xyz.  The  rotation  described  in  Figure  9  on  page  lO.a.  about  the  z  axis  produces 
Figure  10  on  page  10.  The  transformation  from  XYZ  to  xy, Z  is 


x,  =  Ycos  ip  +  F  sin  i p 

y,  =  —X  sin  ip  +  Fcos  ip 
Z  =  Z 


or  written  in  matrix  form 


*i 

Ti 

Z 


COS  ip 

sin  t p 

o' 

'x~ 

—  sin  ip 

COS  i p 

0 

Y 

0 

0 

1 

Z 

(H) 
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Figure  9.  Euler  Angles 


i 

Y 

sin  '/' 

j\ 

COS  1 1> 

COS  <l>  1  'v 

s  \ 

•  sin  i/t 

1  \ 

i 

\ 

Figure  10.  Rotation  about  the  Z  axis  through  the  Angie  \f/ 

Writing  the  transformations  of  Figures  9.b.  and  9.c.  in  matrix  form  yields 


X 

cos  9 

0  -  sin  9 

V 

y\ 

jt. 

sss 

0 

sin  9 

1  0 

0  cost? 

y\ 

_z_ 

(12) 

10 


and 


.T 


y 


1 

0 

0 

X 

0 

cos  4> 

sin  <p 

0 

—  sin  <p 

cos  <p 

_2i_ 

(13) 


Combining  eqs{  11-13)  into  one  equation  yields  the  complete  321  coordinate  transfor¬ 
mation  from  the  inertial  reference  frame  to  the  body  reference  frame 


JC 

~I 

0 

0 

cos  9 

0 

—  sin  9 

COS  ip 

sin  ip 

o' 

y 

= 

0 

cos  <p 

sin  4> 

0 

1 

0 

—  sin  ip 

COS  ip 

0 

J 

_z_ 

0 

—  sin  <p 

cos  <p 

sin  9 

0 

cos  9 

0 

0 

1_ 

K 

cos  9  cos  ip  cos  9  sin  ip 

—  cos  4>  sin  4>  +  sin  <p  sin  6  cos  ip  cos  <p  cos  ip  +  sin  <p  sin  9  sin  ip 
sin  0  sin  i p  +  cos  <p  sin  9  cos  ip  —  sin  0  cos  ip  +  cos  <j>  sin  9  sin  ip 


(14) 
—  sin  9 
sin  </>  cos  9 
cos  <p  cos  9 


X 

Y 

Z 


This  is  just  one  of  the  12  different  Euler  coordinate  transformations  schemes  available. 
Figure  11  on  page  12  lists  all  possible  Euler  angle  combinations  and  their  resulting 
transformation  matrices.  It  should  be  noted  that  eq(  14)  does  not  agree  with  the  321 
Euler  transformation  matrix  of  Figure  11  on  page  12.  The  321  matrix  of  Figure  II  on 
page  12  contains  a  error  that  interchanges  <p  with  \p  and  ip  with  <p . 
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»>rr  .  t 

rin  t  n  M»n  f 
MrntstHiAtiori 


ivrc-i 
ftN.PH  Ann  t 

ncrntstfUAftoM 


c<c# 

-CtMC*  *  1*5* 

sic# 

-itstwctc* 

uuct  ‘  ctu 

u 

-«S# 

etc* 

act 

C4S4C4  *  I >54 

CtUU-1 

ttu 

Jl 

etc* 

ctu 

UC9 

uuct  -c ttt 

ItStSt'CtC* 

CIC# 

SI 

-CIS# 

asm  •  sly 

etc* 

CtUSf 

uct 

S4S#C#  •  CIS# 

-s*e  9 

-S#l#$#  ♦  Cl c# 

etc*  1  uuu 

uct 

-ctu '  uuct 

■uct  •  ctuu 

ceci 

ttu  i  ctuc* 

cist 

-51 

etc* 

act  -5*5454 

C*S4  *  5*54C4 

-uct  ■ 

C454  • 

etc* 

SI 

uct  >  uuu 

im  -ctuu 

C«CI 

cm 

CIS# 

»  1 

-CIS#  •  S*S#C# 

c#c#  •  s#s#s# 

MCt 

5*54  •  C4S4C4 

-t*C4  •  C*5454 

CICI 

5*54 

C*C#  -SICIS# 

ctu  •  ttetet 

C*S# 

-SIC#  -CICIS# 

-UU  •  ctctct 

c# 

SIC# 

UU 

CIS# 

CICIC# -SIS # 

ctctu  •  uct 

1*54 

-side# -cis# 

-nets*  ♦  etc* 

C*C4  -SICIS# 

SIS# 

-ett* -ttetet 

SIS# 

Cl 

MCt 

s*r.4  •  Cf  CIS# 

-CIS# 

-f#S#  •  CICIC# 

CICIC#  -SIS! 

CIS# 

-ClCfS#  -SIC# 

-VC# 

Cl 

Nit 

side#  •  cis# 

SIS# 

-SICIS#  *  c#c# 

C#C#  -SICIS# 

CIS#  •  licit# 

ttu  ■ 

SIC#  -CICIS# 

-sis#  •  tide# 

ctu 

ns# 

-s#c# 

Cl 

CICIC# -SIS# 

CICIS#  •  SIC# 

-ettt 

1 

-S#C#C#  -CIS# 

-f#c#s#  •  act 

MM 

Figure  li.  Table  of  Euler  angle  combinations  )Ref.  4,  p.  764) 
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C.  RIGID  BODY  EQUATIONS 

The  genera]  motion  of  a  particle  and  its  corresponding  momentum  have  been  written 
about  extensively  in  various  publications.  A  summary  of  the  relevent  derivations  and 
equations  of  motion  from  [Ref  6]  are  presented  here. 

1.  Derivative  of  a  Vector  in  a  Rotating  Coordinate  System 

The  time  derivative  of  a  vector  in  a  rotating  coordinate  system  consists  of  the 
rate  of  change  of  the  vector  relative  to  the  moving  axis  and  the  rate  of  change  due  to  the 
rotation  of  the  axis.  Figure  12  shows  a  vector  in  a  rotating  system  where  XYZ  represent 
an  inertial  system  and  xyz  are  the  set  of  axes  rotating  with  angular  velocity  w  relative 
to  XYZ.  If  unit  vectors  along  the  X,  Y,  and  Z  axes  are  i,j,  and  k,  respectively,  the  vector 
r  can  be  written  in  the  form 


r  =  xi  +yj  4-  zk  (15) 


Figure  12.  Vector  in  a  rotating  system  [Ref.  6,  p.  107J 


The  time  derivative  of  r  is 

r  =  xi  +  xi  +yj+yj  +  zk  +  zk  (16) 

From  [Ref.  7],  it  has  been  shown  that  the  time  derivatives  of  the  unit  vectors  are 

/  =  co  x  /  (17) 
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j  =  w  *j 

(18) 

k  =  a)  x  k 

(19) 

Defining  r„,  as  the  time  rate  of  change  of  r  relative  to  xyz,  eq(16)  becomes 

r  =  Kei  +  w  x  r  (20) 


2.  General  Particle  Motion 

Translational  motion  of  a  particle  with  references  to  two  coordinate  systems  is 
shown  in  Figure  13. 


Figure  13.  Translational  motion  of  a  particle  (Ref.  6,  p.  107] 


The  position  vector  of  a  particle  relative  to  coordinate  system  XYZ  is 

R  =  Rc  +  r  (21) 

where  Re  is  the  vector  from  the  axis  of  the  XYZ  coordinate  frame  to  the  xyz  axis  frame 
and  r  is  the  differential  (position)  vector  in  the  xyz  frame.  The  velocity  of  the  particle 
is  the  derivative  of  eq(2l)  with  respect  to  time 

V=R  =  Rc  +  r  (22) 
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The  velocity  at  the  origin  Rc  =  Vc  combined  with  eq(20)  yields 

V  =  Vc  +  Vrel  +  to  x  r  (23) 

3.  Momentum  of  a  Rigid  Body 

A  particle's  momentum  is  the  product  of  its  velocity  and  mass.  The  linear  mo¬ 
mentum  of  the  panicle  of  Figure  13  on  page  14  is 


P=  I  V dm  (24) 

*  m 

where  V  is  the  particle's  velocity  and  m  is  its  mass.  Recalling  that  V„,  is  the  time  rate 
of  change  of  r  relative  to  the  xvz  coordinate  system,  and  assuming  that  the  xyz  axis  is 
centered  in  the  rigid  body,  V„,  =  0.  Substituting  eq(23)  into  eq(24)  with  K,,;=0. 

P  =  I  ( Vc  +  a)  x  r)  dm  =  mVc  +  uy  x  r  dm  (25) 

"m  "m 

By  choosing  the  center  of  mass  of  the  rigid  body  as  the  center  of  the  rotating  coordinate 
system,  f  r  dm  =  0  by  definition.  This  simplifies  eq(25)  to 

P  =  mVc  (26) 

where  P  is  the  linear  momentum  of  the  rigid  body,  m  is  its  mass,  and  Vc  is  the  velocity 
of  its  center  of  mass. 

The  angular  momentum  of  a  particle  dm  about  point  C  is  defined  as  the  moment 
about  C  of  the  linear  momentum  given  by 

hc  =  r  x  V  dm  (27) 


Substituting  for  V  yields 


Hc  =  I  r  x  V  dm  =  r  x  (Vc  +  co  x  r)  dm 

•'m 


Rearranging  terms 


(28) 
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Hc  =  |  f  r  dm  \  x  Vc  +  j  rx(wxr)  dm  (29) 

X*'  m  J  "  m 

By  choosing  the  center  of  mass  of  the  rigid  body  as  the  origin  of  the  rotating  coordinate 
system,  eq(29)  reduces  to 


Hc  =  I  r  x  (co  x 

*  m 


r)  dm 


(30) 


Reducing  the  angular  velocity  vector  co  into  its  components 


co  =  (oxi  +  coy  +  co.k 


(3i: 


Combining  equations(15),  (30),  and  (31)  yields 


Hc  =  [Jxx^x  ~  y,(t)y  —  ^X2W2^1 

+  [  lxyCJ  x  +  lyytOy  l y  jlO  d\) 

+  [  —IXMX  —  I yyCOy  +7^.(0^]^ 


(32) 


where 


J  m 

=  J  (. x2  +  z 2) 

J  m 

-J  (^2+/) 


dm 


dm 


dm 


are  the  moments  of  inertia  about  the  x,y,  and  z  axes,  respectively,  and 


4y  =  J  xydm 


(33) 


(34) 


(35) 


(36) 
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=1 


xz  dm 


yz  dm 


(37) 


(38) 


are  the  products  of  inertia.  Equation(32)  written  in  matrix  form  becomes 


Hx 

Hy 

= 

H, 

7  _  r 

lxx  lxy 


-1 


f 

XZ 


lyy 


-4 


y: 


OJr 


0)v 


(39) 


4.  Equations  of  motion/ Euler's  Moment  Equations 

The  derivation  of  the  following  Euler's  moment  equations  are  summarized  from 
(Ref  6].  Recalling  that  angular  momentum  hc  relative  to  the  point  C  is 


hc  =  r  x  ntv 


(40) 


its  derivative  with  respect  to  time  is 

hc  =  r  x  mv  +  r  x  mi  (41 ) 

The  moment  of  the  force  about  a  point  C  by  definition  is 

Me  =  r  x  F  (42) 

Recall  Newton's  second  law  describing  force: 

F  =  ma  =  mv  (43) 

Substituting  eqs(22)  and  (43)  into  eq(41)  yields 

hc  =  rx  m(Rc  +  r)  +  r  x  F  (44) 

Substituting  eq(42)  into  eq(44)  and  rearranging  terms  yields  the  moment  equation 


17 


MC  =  HC  + 


( Rc  x  r)  dm 


(-45) 


where  Hc  is  the  angular  momentum  of  the  rigid  body  about  point  C,  From  eq(20)  the 
latter  part  of  eq(45)  can  be  written  as 


Rcx  r  dm  =  Rc  x 

*  m 


(46) 


Noting  that  point  C  is  the  center  of  mass  f  r  dm  =  0  and  r„,  =  0  for  a  rigid  body,  eq(46) 
reduces  to 


Mc  =  Hc  (47) 

The  moment  of  external  forces  about  the  center  of  mass  of  a  rigid  body  is  equal  to  the 
time  rate  of  change  of  the  angular  momentum  of  the  body  about  the  center  of  mass  [Ref. 
7|.  In  terms  of  a  rotating  axis,  the  moment  equation  of  eq(47)  becomes 


Mc  —  Hcrel  +  co  x  Hc  (4$) 

Substituting  eq( 39)  in  for  H.  the  components  of  the  moment  equation  can  be  written  as 

Mx  —  Hx  +  c oyHz  —  u)2Hy  (49) 

My  =  Hv  +  <o2Hx  —  (oxH t  (50) 

Mz  =  H2  +  ( oxHy  —  (oyHx  (51) 


where  are  the  moment  equations.  In  this  study,  an  assumption  is  made  that  the 
principle  axes  of  the  moment  of  inertia  are  the  rotating  axes  x,y,z;  therefore,  the  pro¬ 
ducts  of  inertia  (/„,  /„,  Iyt)  are  equal  to  zero.  The  angular  momentum  components  re¬ 
duce  to 

Hx  =  Ixx^x  (52) 
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Hy  lyy(Oy 


(53) 


H2 


(54) 


Substituting  eqs(52-54)  into  the  moment  equations  yields  Euler's  moment  equations 

MX  ~  ^XX^X  “b  t0yC02(  I 22  lyy)  (55) 

yiy  =  I yylOy  4"  CO xCO 2i  I XX  I 22^  (56) 

M:  =  l22U>2  +  "jc <°yUyy  ~  4x)  (57) 


Referring  to  the  coordinate  system  shown  in  Figure  6  on  page  5,  the  angular 
velocities  in  terms  of  orbital  rate  co0  and  attitude  error  angles  pitch(0),  roll(0),  and  yaw 
(0)  from  [Ref.  6,  p.  130]  are 


CO* 

1 

0 

—  sin  6 

0 

C Oy 

= 

0 

COS  0 

cos  d  sin  0 

e 

0 

—  sin  0 

cos  6  cos  0 

0 

cos  6  sin  i Jj 

cos  0  cos  0  +  sin  0  sin  6  sin  0 
—  sin  0  cos  0  4-  cos  0  sin  6  sin  0 


(5S) 


When  using  small  angle  approximations,  eq{58)  reduces  to 


CO* 

0  ~  COO0 

(Oy 

= 

6 -co0 

0  +  COO0 

Note  that  the  orbital  rate  is  defined  as 


(59) 


(60) 


where  n  is  the  gravitational  constant  and  R  is  the  radius  of  the  earth  plus  the  satellite's 
altitude.  For  the  purpose  of  station  keeping,  the  small  angle  approximations  are  valid. 
Substituting  eq(59)  into  the  rigid  body  moment  eqs(55-57)  yields  the  linearized  rigid 
body  dynamic  equations 
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(6!) 


Mx RB  —  ^xx&  +  <*>o(Iyy  ~  —  a)o(^«  +  hz  ~  lyy)^ 

MyRB  =  lyye  (62) 

^zRB  =  I 22^  ^oUyy  ^xx)^  ^O^xx  "b  ^zz  Jyy)*P  (63) 


D.  GRAVITY  GRADIENT  EFFECT 

Attitude  control  using  gravity  gradient  stabilization  works  on  the  principle  that 
given  an  object  with  an  asymmetrically  distributed  mass,  the  object  will  tend  to  align  it¬ 
self  with  its  minimum  moment  of  inertia  axis  along  the  local  vertical  towords  the  earth. 
The  force  due  to  gravity  on  a  mass  element  is 


dF  =  - 


R 3 


R 


(64) 


w'here 


and 


R  =  Rq  +  T 


He  =  =  SeRe 


(65) 


(66) 


where  G  =  universal  gravity  constant,  M,~  the  mass  of  the  earth,  g,=  gravity  acceler¬ 
ation  at  the  earth's  surface,  and  R,  =  the  radius  of  the  earth.  Reference  8  shows  the 
torque  generated  about  the  mass  center  of  the  object  is 


d\lg  =  r  x  dF  = 


Hfdm 

R3 

Hedm 

R 3 

Hgdm 

R3 


r  x  R 


r  x  (Rq  +  r) 


(67) 


rxRo 


Cubing  eq(65)  yields 


R '  =  Rn 


,  ,  3'«o  ,  ,  r3 

1  +  >2  '  n2  +  r>3 


*0 


*0 


(68) 
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Noting  that  r<$Ra,  the  higher  order  terms  drop  out.  Inversion  of  eq(68)  yields 


11/.  «  Rp 

R3  Ro3  V 


(69) 


Substituting  eq(69)  into  eq(67)  and  integrating  to  find  the  moment  yields 


M. 


=  -  He  [  r  dm  x  Rq  +  [ 

Rq 


f  (r  x  Ro)(r .  Rq) 


dm 


(70) 


Given  that  the  center  of  mass  is  the  origin  of  the  body  frame,  f  r  dm  =  0,  eq(70)  reduces 
to 


A/,— (rxRo)(''Ro)dm  (71) 

R0 


In  the  reference  system  defined  by  Figure  14,  the  K  (Z)  vector  is  pointing  towards 
the  earth;  therefore,  the  vector  R„  in  the  reference  coordinate  frame  is 

*or«*oK  (?2) 
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Restating  the  Euler  321  rotation  matrix  described  earlier 


/ 

’l 

0 

0 

cos  6 

0 

—  sin  6 

COS  0 

sin  0 
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j 
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0 

COS  0 
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0 
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0 

—  sin  0 

COS  0 

0 

j 

_k~ 

0 

—  sin  0 

COS  0 

sin  6 

0 

cos  6 

0 

0 

1_ 

K 

cos  0  cos  0  cos  0  sin  0 

—  cos  0  sin  0  4-  sin  0  sin  0  cos  0  cos  0  cos  0  +  sin  0  sin  0  sin  0 
sin  0  sin  0  +  cos  0  sin  0  cos  0  —  sin  0  cos  \p  +  cos  <p  sin  0  sin  i p 


(73) 

-i 


—  sin  6 

sin  0  cos  6 

J 

cos  0  cos  6 

K_ 

Transforming  R0,  into  the  spacecraft  s  fixed  coordinates  i,  j,  and  k: 


cos  9  cos  0 

cos  0  sin  0 

—  sin  0 

‘o' 

7?0  r  ~ 

—  cos  0  sin  0  +  sin  0  sin  0  cos  0 

cos  0  cos  0  +  sin  0  sin  0  sin  0 

sin  0  cos  0 

0 

sin  0  sin  0  +  cos  0  sin  0  cos  0 

—  sin  0  cos  0  +  cos  0  sin  0  sin  0 

COS  0  cos  0 

A' 

which  reduces  to 


R0 r  =  R0{  -  sin  0/  +  sin  <p  cos  6j  +  cos  <p  cos  6k) 
Substituting  eq(75)  into  eq(71)  yields  the  gravity  gradient  moment  equations 

{I22  —  />v)  sin  <p  cos  <p  cos L6 

—  —  /7Z)  sin  6  cos  6  cos  <p 

—  ( Iyy  —  lx*)  sin  6  cos  6  sin  </> 


Replacing 


(75) 


(76) 


(77) 


where  co0  is  the  orbital  rate,  and  making  the  small  angle  approximations,  the  gravity 
gradient  moment  equations  become 


{J 22  lyy)4* 

(/„  -  U0 


0 


(78) 
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£quation(78)  clearly  shows  that  by  increasing  the  moments  of  inertia  for  7,  and  It 
with  respect  to  /,,  the  gravity  gradient  moments  for  roll  and  pitch  will  increase.  By  ex¬ 
tending  a  point  mass  along  the  z  axis,  a  sizeable  increase  in  the  moments  of  inertia  along 
the  x  and  y  axes  is  realized. 


The  magnitudes  of  the  gravity  gradient  torques  generated  at  different  altitudes  is 
given  by  Figure  15. 


Figure  15-  Gravity  Gradient  torque  Magnitudes  vs  Altitude  [Ref.  3,  p.  58] 

E.  DESIGN  FOR  STABILITY 

1.  Stability  and  Frequency  Analysis 

Combining  the  rigid  body  dynamic  equations  (61-63)  with  the  gravity  gradient 
equation  (78)  yields  the  undamped  linear  dynamic  equations  for  roll,  pitch,  and  yaw: 

Roll  lj>  +  Acol(Iy  -  I2)4>  -  c o0{Ix  +  1,-  Iy)4 '  =  0  (79) 

Pilch  ly8  +  3(00(7*  -  I,)0  =  0  (80) 
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(81) 


Yaw  I2\p  +  cosily  -  Ix)\jj  +  co0{Ix  +  4  —  Iy)4>  =  0 
Using  Reference  9  notations  for  simplicity,  the  following  variables  are  defined: 


(82) 

(83) 


(84) 


Substituting  eqs( 82-84)  into  eqs( 79-81)  yields  the  following  second  order  differential 
equations 


< i i  4-  AA(ol<j>  —  (1  —  A)co0i[/  =  0 

(85) 

£  +  lBco20d  =  0 

(86) 

+  Ccolip  +  (1  —  Qc(>  =  0 

(87) 

Taking  the  Laplace  transform  of  eqs(85-87)  allows  for  frequency  analysis.  The  results 
are  put  into  matrix  form: 


S2  +  4  A  col  0 

—  (1  —  A)co0S 

V 

‘o’ 

0  S2  +  3Ba>l 

0 

e 

= 

0 

(1  -  OcooS  0 

S2  +  Ccol 

A 

0_ 

(88) 


The  determinant  of  the  matrix  is  the  system's  characteristic  equation.  Using  the  middle 
row  to  simplify  calculations,  the  characteristic  equation  is  found  to  be 

(S4  +  ( 1  +  3A  +  CA)co20S2  +  4  A  C«o)(S2  +  35wo)  =  0  (89) 

“Solution  of  the  first  term  on  the  left  yields  the  natural  frequencies  of  roll  and  yaw, 
repectively;  the  second  term  on  the  left  defines  the  natural  frequency  of  the  pitch  axis." 
(Ref.  9,  p.  9]  Solving  for  the  first  term  on  the  left  (roll, yaw  stability)  by  the  quadratic 
formula  and  dividing  through  by  col  yields 
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(90) 


S20l,  Sf,aw  -(l  +  3A  +  AC)±^(l  +  3A  +  AC)2-\6AC 

2  ~  2  ~  ■> 

fc>o  co0 

For  roll  yaw  stability  the  roots  of  eq(90)  must  be  real  and  negative;  therefore  the  fol¬ 
lowing  stability  conditions  exist  [Ref.  4  p.  611] 

l  +  3A  +  AC>4j~AC  (91) 


AC>  0 


Solving  for  the  second  term  of  Eq(89)  and  dividing  through  by  yields 


=  35  = 


3(C  —  A) 
1  -  CA 


(92) 


(93) 


The  condition  for  stability  for  the  pitch  axis  is 

A  >  C  (94) 

Converting  A.  B.  and  C  back  into  their  respective  moment  of  inertia  relationships  de¬ 
fined  by  eqs(  82-84).  combined  with  the  conditions  for  stability,  yields  two  possible  ori¬ 
entations:  /,>/,>  /z  and  /,>/,>  Iy.  Figure  16  on  page  26  shows  a  summary  of  the 
stability  regions  based  on  moment  of  inertia  ratios.  [Ref.  10  p.  204] 
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Figure  16.  Gravity  Gradient  Stability  Regions 
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III.  ENVIRONMENTAL  DISTURBANCES  AND  RESTORING  TORQUES 

A.  ENVIRONMENTAL  DISTURBANCES/TORQUES 

The  are  four  major  environmental  disturbances  that  affect  an  orbiting  spacecraft: 
aerodynamic,  magnetic,  gravity  gradient,  and  solar  radiation  torque.  Figure  17  lists 
these  disturbances  and  their  relative  regions  of  dominance. 
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Figure  17.  Environmental  Disturbance  Torques  [Ref.  4,  p.  17] 


1.  Aerodynamic  disturbance 
a.  Drag 

“For  satellites  near  the  earth  the  principle  non-gravitational  force  is  aero¬ 
dynamic  drag.  Aerodynamic  drag  is  a  retarding  force  due  to  atmospheric  friction  and  is 
in  the  direction  opposite  the  space  vehicles  velocity  vector."  [Ref.  4,  p.  63] 

The  drag  force  is  a  function  of  vehicle  velocity,  air  density,  and  surface  area 
of  the  satellite.  The  drag  force  is  given  by  the  following  equation: 


2 


(95) 


where  d  is  the  drag  force,  p  is  the  atmospheric  density,  v  is  the  velocity  of  the  satellite 
along  its  orbital  path,  c4  is  the  drag  coefficient,  and  a  is  the  area  of  the  satellite  over 
which  the  drag  acts. 

b.  Aerodynamic  torque 

The  collisions  of  air  molecules  of  the  upper  atmosphere  with  the  satellite 
surface  produces  a  torque  about  its  center  of  mass.  Aerodynamic  disturbances  are  a 
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function  of  the  satellite  s  altitude,  velocity,  and  symmetry.  A  satellite  with  a  relatively 
low  earth  orbit  may  be  significantly  affected  by  aerodynamic  torques. 

Shown  in  Figure  18  is  a  cylindrical  satellite  with  its  center  of  mass  (cm)  and 
center  of  pressure  (cp).  The  offset  between  the  center  of  mass  (cm)  and  the  center  of 
pressure  (cp)  is  denoted  by  L„  Lf,  and  L„  respectively.  The  simplified  aerodynamic  dis¬ 
turbance  torques  from  Ref  8  may  be  computed  by  the  following  expression: 

Ta  =  pJaAa  sin  a  (96) 

where  T ,  is  the  torque  due  to  aerodynamic  pressure  (ft-lbs),  P,  is  the  aerodynamic  pres¬ 
sure  (lbs,'//J),  4  is  the  distance  between  the  center  of  mass  and  the  center  of  pressure  (ft), 
A „  is  the  exposed  surface  area  (ft2)  and  a  is  the  angle  of  attack  (radians).  For  practical 
purposes,  the  aerodynamic  torque  in  the  x  direction  is  zero  (sin  a  =  0  when  ocsO).  The 
conversion  factor  from  ft-lbs  to  N-m  is  approximately  1.356.  Figure  19  on  page  29 
shows  aerodynamic  pressure  as  a  function  of  altitude. 


Figure  18.  Aerodynamic  effects  for  cylindrical  satellite  [Ref.  3,  p.  58] 
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Figure  19.  Aerodynamic  pressure  as  a  function  of  altitude  [Ref.  8,  p.  456] 


The  magnitude  of  the  pitch  and  yaw  aerodynamic  disturbances  for  a  one 
inch  cm-cp  offset  are  shown  in  Figure  20  on  page  30  for  various  altitudes.  The  magni¬ 
tude  of  these  torques  decrease  significantly  with  increases  in  altitude.  "Proper  control  of 
the  satellites  inertias  to  enhance  gravity  gradient  torques  can  make  aerodynamic  dis¬ 
turbances  essentially  insignificant  for  orbital  altitudes  above  300  nautical  miles  (480km). 
On  the  other  hand,  in  the  100  to  200  nautical  mile  (160-320km)  altitudes  aerodynamic 
disturbances  are  3  to  4  orders  of  magnitude  greater  than  the  gravity  gradient  torque  from 
a  1  slug-fr1  inertia  difference  and  1  degree  attitude  error.  In  the  lower  regions,  very  large 
inertia  differences  are  required  for  accuracy.  [This  leads  to  design  considerations]  where 
proper  configuration  can  make  these  aerodynamic  torques  work  as  restoring  torques 
instead  of  disturbance  torques"  [Ref.  3,  p.  57]. 
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Figure  20.  Aerodynamic  torque  for  a  1  inch  offset  as  a  function  of  altitude. 

2.  Rigid  Body  Solar  Pressure  Torque 

"The  major  factors  determining  the  radiation  torque  on  a  spacecraft  are  the  in¬ 
tensity  and  spectral  distribution  of  the  incident  radiation,  the  geometry  of  the  surface 
and  its  optical  properties,  and  the  orientation  of  the  sun  vector  relative  to  the  spacecraft" 
[Ref.  4,  p.  570]. 

Solar  pressure  torques  are  the  disturbances  produced  by  solar  radiation  pres¬ 
sure  and  is  a  function  of  the  offset  of  the  center  of  pressure  (cp)  and  the  center  of  mass 
(cm).  Referring  to  Figure  21  on  page  31,  the  center  of  mass-center  of  pressure  offsets 
are  designated  as  L,,  Ly,  and  L„  respectively.  The  torques  generated  due  to  solar  pres¬ 
sure  are  given  by  (Ref.  8,  p.  454J  as: 


Ts  =  PSLSAS  cos  e  (97) 

Ts  =  2 PSLSAS  cos2e  (98) 

where  T,  =  torque  due  to  solar  radiation  (ft-lbs),  P,  -  radiation  pressure  (lb/fir2)  and  is 
»  constant  at  9.65  x  10'*  lb !fP  for  an  earth  orbiting  vehicle,  L,  =  center  of  mass  center 
of  pressure  offset  (ft),  A,  =  surface  area  of  satellite  normal  to  sun(/P),  and  e  =  angle 
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of  incidence  (degrees).  Note,  equation  (97)  is  for  an  absorbent  satellite  body  and 
equation  (98)  is  for  a  reflective  body. 


V«-\ 

vc- 

M 

n 

i 

i 

,  A 

L, 

— 

1 

’ 

Figure  21.  Rigid  body  solar  pressure  torques  [Ref.  3,  p.  61J 


B.  RESTORING  TORQUES 

As  seen  from  chapter  2,  the  restoring  torques  of  the  undamped  rigid  body  combined 
with  the  gravity  gradient  effects  are  considerably  smaller  for  the  yaw  axis  than  for  either 
the  pitch  or  roll  axes.  This  is  due  to  the  necessity  of  large  moments  of  inertia  along  the 
pitch  and  roll  axes,  as  compared  to  the  yaw  axis,  for  gravity  gradient  stabilization  to 
occur.  The  lack  of  adequate  yaw  restoring  torque  requires  that  additional  yaw  restoring 
torque  be  provided  for  accurate  3-axis  stabilization.  A  discussion  of  yaw  restoring  torque 
methods  follows. 

1.  Thrusters 

Thrusters  are  an  effective  torque  generating  device;  however,  for  the  purpose 
of  this  thesis,  thrusters  will  not  be  considered  for  attitude  control.  Thrusters  are  neces¬ 
sary  for  orbit  maintenance  and  orbital  transfer.  Due  to  the  physical  size  constraints  im¬ 
posed  by  the  GAS  cannister,  ORION  is  severely  limited  in  its  fuel  storage  capacity  to 
approximately  71  lbs  of  hydrazine  fuel  [Ref  1J.  An  assumption  that  the  launch  platform 
will  place  ORION  in  its  final  orbit  is  optimistic.  Figure  22  on  page  32  shows  the  fuel 
required  to  go  from  one  orbital  altitude  to  another.  Saving  the  onboard  fuel  for  orbit 
maintenance  and  transfer  and  not  for  attitude  control  allows  more  fuel  to  be  devoted 
towards  supporting  the  mission  either  by  allowing  greater  orbital  flexibility  and/or  a 
longer  lifetime. 
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2.  Momentum  Exchange  Devices 

The  principle  behind  momentum  exchange  devices  is  that  by  changing  the  an¬ 
gular  momentum  of  the  control  device,  the  angular  momentum  of  the  vehicle  will  be 
changed  an  equal  and  opposite  amount.  Recall  that  the  rate  of  change  of  angular  mo¬ 
mentum  is  equal  to  the  torque  that  is  generated.  The  two  basic  types  of  momentum  ex¬ 
change  devices  are  reaction  wheels  and  control  moment  gyros  (CMG).  The  basic 
difference  between  the  two  is  the  method  by  which  the  change  in  angular  momentum  is 
accomplished.  The  reaction  wheel  axis  is  fixed  (generally  along  one  of  the  axes  of  the 
vehicle)  so  that  the  change  in  angular  momentum  is  accomplished  by  varying  the  speed 
of  rotation  of  the  flywheel.  The  CMG  creates  the  change  in  angular  momentum  by  tilt¬ 
ing  the  contant  speed  flywheel  with  respect  to  the  vehicle. 
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a.  Reaction  wheels 


"An  axis  may  be  controlled  by  varying  the  speed  of  the  reaction  wheel  in 
response  to  an  attitude  error  [Ref.  6,  p.  149].  The  momentum  of  a  system  with  a  reaction 
wheel  is 


H=HW  +  HC 


where  Hc  is  the  angular  momentum  of  the  system  and 


Hw  = 


Hwx 

Hwy 

Hw: 


(99) 


is  the  angluar  momentum  of  the  reaction  wheel.  Aligning  the  reaction  wheels  along  each 
axis,  and  substituting  eq(99)  into  eq(4S)  and  adding  the  gravity  gradient  moments  of 
eq(7S)  yields 


40  4-  4a»o(fy  ~  4)0  ~  &>o(4  +  4  “  4)0  +  ~  -  (0  +  "o 4>)Hwy  +  Hwx  =  X!7* 

Iy6  +  3 Wq(/x  -  4)0  +  (  -  4>  +  o0\]j)Hw2+  (iA  4-  u)0(f))Hwx  4-  Hwy  -  X Ty  (101) 

40  +  "o(4  -  4)0  +  wo(4  +  4  -  4)0  +  (0  -  ^o0)  +  ( -  0  +  wo0)J4*  +  Hw2  =  £T2 

"Because  Hw„HwrH„,  and  <y0  are  small,  the  coupling  terms  are  small.  If  the  coupling 
terms  are  neglected,  the  equations  of  motion  about  the  roll,  pitch,  and  yaw  axes  become 
independent,  and  hence,  they  can  be  controlled  independently"  (Ref.  6,  pp.  149-150], 
Dropping  the  coupled  terms,  the  motion  equations  become 

40  4-  4a»o(4  ~  4)0  “  wo(4  +  4  “  4 )0  +  HWx  =  X7*  (,02) 


iyd  +  3<»o(4  —  4)0  +  =  X4 


(103) 


40  +  wo(4  ~  4)0  +  wo(4  +  4  -  4)0  + 


(104) 
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Allowing  the  momentum  of  the  reaction  wheel  to  be  proportional  to  the  attitude  error, 
then  by  differentiation,  the  reaction  wheel  control  torques  are  equal  to 

T„w  =  ~  ^(^0  +  0)  (105) 


Tyrw  =  “  +  0) 


(106) 


7,™  =  -  A'2(r7^  +  (107) 


and  the  control  scheme  would  be  as  shown  in  Figure  23. 


Figure  23.  3  axis  reaction  wheel  control  scheme  |Ref.  6,  p.  149] 


The  torque  genterated  by  the  reaction  wheel  is  a  function  of  the  wheel's  angular  velocity 
and  its  moment  of  inertia.  This  means  that  a  small  lightweight  wheel  at  a  high  rate  of 
speed  generates  the  same  amount  of  torque  that  a  large  wheel  at  a  low  rate  of  speed. 
Since  ORION  is  designed  for  a  relatively  short  lifespan,  wear  and  tear  on  spinning  parts 
is  not  an  overriding  concern.  It  should  be  noted  that  continuous  disturbance  torques 
that  are  not  cyclical  in  nature  will  eventually  saturate  the  reaction  wheel,  whereby,  the 
reaction  wheel  rotational  speed  will  reach  a  maximum.  Thrusters  or  some  other  method 
of  momentum  removal  must  be  used  to  reduce  wheel  speed.  Disturbance  torques  that 


34 


are  cyclical  in  nature  are  easily  olTset  by  the  reaction  wheel  as  long  as  the  disturbance 
torque  is  such  that  it  does  not  saturate  the  reaction  wheel  speed. 

3.  CMG(s) 

"A  control  moment  gyro's  (either  single  or  double  gimbaled)  angular  mo¬ 
mentum  is  due  to  the  rotor  which  is  spinning  about  the  spin  axis  with  a  constant 
angular  rate.  Because  the  spin  axis  is  gimbaled,  a  commanded  gimbal  rotation  causes 
the  direction  of  the  angular  momentum  vector  to  change,  thus  creating  a  control 
torque  parallel  to  the  output  axis.  The  magnitude  of  this  torque  depends  on  the 
speed  of  the  rotor  and  the  gimbal  rotation  rate."  [Ref.  4,  p.  200] 


A  control  moment  gyro  scheme  is  shown  in  Figure  24. 


Figure  24.  Twin  gyro  controller  with  one  drive  motor  [Ref.  8,  p.  415] 


Because  of  their  expense  and  weight,  CMG's  are  only  used  on  large  spacecraft 
[Ref.  4],  For  this  reason,  control  moment  gyros  will  not  be  addressed. 

4.  Magnetic  torquers 

Magnetic  torquers  may  be  used  to  dampen  the  kinetic  energy  of  the  libration 
motion  as  well  as  correcting  attitude  errors  in  any  direction.  The  magnetic  torquers 
consist  of  three  wire  coil  sets  arranged  perpendicularly  to  each  other.  This  configuration 
allows  magnetic  torque  to  be  generated  along  any  axis  or  in  any  direction  as  a  result  of 
adding  two  or  three  different  orthogonal  torques.  As  shown  in  Figure  25  on  page  36, 
w  hen  a  current  is  applied  around  the  loop  of  a  coil,  it  produces  a  magnetic  dipole  that 
is  normal  to  the  plane  of  the  coil  and  a  magnitude  which  is  proportional  to  the  coil's 
enclosed  area  and  ampere-turns. 


Figure  25.  Magnetic  Moment  Due  to  a  Current  Loop  [Ref.  4,  p.  204] 


I  he  magnetic  moment  M  is  defined  as 

M  =  XI  An  (108) 

where  N  is  the  the  number  of  turns  around  the  coil,  I  is  the  current  in  amperes,  A  is  the 
area  enclosed  by  the  loop,  and  h  is  the  unit  vector  normal  to  the  plane  of  the  loop.  Be¬ 
cause  of  the  relationship  between  N,  1,  and  A,  many  different  combinations  are  available 
to  produce  the  desired  magnitude  of  the  magnetic  moment.  Figure  26  on  page  37  shows 
three  existing  spacecraft  with  various  combinations  of  number  of  turns  and  current 
drawn.  The  torque  generated  by  the  magnetic  torquer  is 

Tmag  =  MxB  (109) 

where  M  is  the  magnetic  dipole  and  B  is  the  earth's  magnetic  field.  The  earth's  magnetic 
field  for  an  equatorial  orbit  can  be  approximately  modeled  as  a  tilted  dipole  located  at 
78.3°  north  latitude  and  69°  west  longitude  as  shown  in  Figure  27  on  page  37. 
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SPACECRAFT 


REMARKS 


SAS-3  COIL  CONSISTS  OF  260  TURNS.  MAXIMUM 

CURRENT  IS  0.6  A  AND  MAXIMUM  POWER 
CONSUMPTION  IS  10  W 

OSO-8  COIL  CONSISTS  OF  360  TURNS.  MAXIMUM 

CURRENT  IS  0.075  A 

AE-3  TWO  COILS;  EACH  HAS  500  TURNS.  MAXI¬ 

MUM  POWER  CONSUMPTION  IS  12  W 


Figure  26.  Some  existing  magnetic  moment  schemes  [Ref.  3,  p.  205] 


Figure  27.  Earth's  magnetic  Field  modeled  as  a  tilted  dipole  [Ref.  6,  p.  148] 


The  strength  of  the  earth's  magnetic  field  Tor  a  circular  equatorial  orbit  in  component 
form  can  be  approximated  by 
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(110) 


M. 

Bx  - - j-  sin(l  1.7)  sin(281  —  lams) 

R 

Mg 

By  = - f-  cos(l  1.7)  (111) 

R 


Mg 

B2  = — f-  sin(l  1.7)  cos(281  —  /am,)  (112) 

R 

where  lams  is  measured  in  degrees  of  right  ascension  from  Grenich,  Mt  is  the  strength 
at  the  earth's  center  which  is  approximately  8  x  I01S  Wb-m,  and  R  is  the  distance  in 
meters  from  the  earth's  center  to  the  orbital  altitude.  With  respect  to  the  satellite,  the 
earth's  magnetic  field  is  changing  continuously  as  the  vehicle  moves  along  its  orbital 
path.  Putting  into  component  form  both  M  and  B,  eq(109)  becomes 


Tmag-x  ~  M2By  —  MyB2 

(113) 

Tmag-y  =  MxB2  ~  M2Bx 

(114) 

Tmag-2  =  S1yBx  —  MxBy 

(115) 

It  should  be  noted  that  the  earth's  magnetic  field  is  not  uniform  and  varies  in 
direction  based  on  the  satellite's  position.  For  accurate  attitude  control,  the  vehicle  must 
be  able  to  exactly  know  both  its  position  and  the  earth's  magnetic  field  at  that  position. 
For  low  earth  orbits,  the  magnetic  field  is  well  mapped  and  may  be  easily  programmed 
into  the  attitude  control  computer.  For  the  purpose  of  this  thesis,  magnetic  field  ap¬ 
proximations  will  suffice  to  demonstrate  magnetic  torquing  as  a  viable  attitude  control 
capability. 
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IV.  CHAPTER  4 


A.  GRAVITY  GRADIENT  STABILIZATION  RESULTS 
1.  Attitude  Control 

There  are  three  general  attitude  control  phases  with  respect  to  gravity  gradient 
systems:  gravity  gradient  capture,  transition  period,  and  steady  state.  Gravity  gradient 
from  a  tumbling  mode  occurs  when  angular  velocity  slows  to  one  rotation  per  orbit  or 
less,  allowing  the  minimum  moment  of  inertia  axis  to  align  with  the  local  vertical.  The 
transition  period  is  the  time  from  gravity  gradient  capture  until  a  steady  state  has  been 
achieved.  This  period  of  oscillations  will  vary  depending  on  the  dampening  scheme  used 
to  reduce  the  amplitude  of  oscillations.  Predominant  throughout  the  satellite's  lifetime, 
the  steady  state  phase  is  the  period  in  which  the  satellite's  errors  are  confined  within 
operationally  acceptable  amplitudes.  Maintenance  of  the  desired  steady  state  requires 
minimization  of  external  disturbances  and  control  of  those  which  cannot  be  eliminated. 

The  following  group  of  simulations  assumes  that  a  steady  state  phase  has  been 
achieved  and  that  the  small  angle  approximations  are  in  effect.  For  each  of  the  simu¬ 
lations,  the  following  initial  conditions  are  in  effect:  4>,  9,  and  i/>  position  error  =  5°,  and 
</>,  6  ,  and  <//  =  0. 


2.  Various  results  based  on  different  configurations 

a.  Effects  of  extending  booms 

For  Figure  28  on  page  40,  a  10  kg  point  mass  is  used  and  the  mass  of  the 
boom  is  assumed  to  be  negligible.  Extension  of  the  boom  increases  the  moment  of  in¬ 
ertia  in  the  x  and  y  directions.  This  increases  the  restoring  torque  along  the  roll  and  pitch 
axes  although  it  does  little  to  dampen  the  libration  motion. 

b.  Gravity  gradient  effect  alone 

The  following  figures  show  the  gravity  gradient  effect  on  a  satellite  at  an 
altitude  of  1000  km  having  moments  of  inertias  as  shown.  As  shown  by  Figure  29  on 
page  41  and  Figure  30  on  page  42,  stabilization  along  the  roll  and  pitch  axes  would  take 
a  prohibitively  long  time  while  the  yaw  direction  would  not  stabilize.  It  is  clear  that 
gravity  gradient  stabilization  alone  is  incapable  of  acceptable  three  axis  stabilization. 


MOMENT  OF  INERTIA 


MOMENT  OF  INERTIA  VS  DOOM  LENGTH 
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Figure  30.  Gravity  gradient  effect  only  Ix=  95,  Iy=  98,  lz=  3 


c.  EJJccts  of  Reaction  \Wteel(  s) 

One  method  of  augmenting  gravity  gradient  stabilization  is  through  the  use 
of  reaction  wheels.  Teldrix  reaction  wheel  classes  and  corresponding  technical  data  are 
presented  in  Figure  31  on  page  43. 
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Wheel  diameter 

cm 

20 

26 

35 

50 

60') 

Angular  momentum  range 

Nms 

1  0  ...6.5 

5.0...  20 

14. ..00 

50... 300 

200...  1000 

Max.  ronclion  torque 

Nlll 

0.2 

0.2 

0.2 

0.3 

0.3-0. 6 

Speed'") 

min  1 

6000 

6000 

6000 

6000 

6000 

Loss  torqrrc  at  mart.  speed") 
Power  consumption: 

Nm 

<0.012 

<  0.013 

<  0.015 

S  0.022 

<  0.07 

-  steady  stale  (depending 
on  speed) 

w 

2. ..7 

2...B 

2... 10 

3... 16 

10. .50 

—  max.  power  rating 

w 

<  60 

<  00 

<  100 

<  150 

<  500 

Dimensions: 

-  diameter  A 

mm 

203 

260 

350 

500 

600 

-  height  B 

mm 

75 

05 

120 

150 

180 

Wfcight 

2. 7.. .3.4 

3.5.. .6.0 

5.0... 0.0 

7.5... 12 

20. .37 

Environmental  conditions: 

-  operating  temperature 

-  vibration  (sinusoidal) 
-vibration  (random) 

-  linear  acceleration 

suitable  lor  satellites  compatible  with  Inunchc'S 
such  as  ARIANE  or  Space  SI  rutile 

Figure  31.  Teldrix  Reaction  Wheel  Technical  Data 


Of  the  five  diameter  classes,  the  20  cm  diameter  wheel  is  the  most  appropriate  of  the 
classes  listed  for  ORION.  Reaction  wheels  also  come  in  diameters  smaller  than  20  cm 
diameter.  Bendix  Aerospace  SP  reaction  wheel  of  Figure  32  on  page  44  has  a  diameter 
of  6  inches. 

As  illustrated  in  Figure  31  and  Figure  32  on  page  44,  power  consumption 
values  differ  with  dilferent  reaction  wheels.  For  the  20  cm  Teldrix  reaction  wheel,  steady 
state  requires  2-7  watts  with  a  maximum  less  than  60  watts.  The  Bendix  SP  reaction 
wheel  requires  2.5  watts  at  steady  state  with  a  maximum  power  consumption  of  6.5 
watts.  Due  to  weight  and  power  consumption  values,  the  Bendix  reaction  wheel  appears 
to  be  the  superior  choice  of  the  two. 

The  following  figures  were  generated  by  setting  the  reaction  wheel  torques 
equal  to  eqs(107-109).  The  initial  conditions  are  <f>  ,  8,  and  ^  position  errors  =  5°,  and 
(f)  =  6  =  t/r  =  0.  For  simulation  purposes,  the  maximum  reaction  torque  value  of  .2  N-m 
from  Figure  31  will  be  in  effect. 

Figures  33-35  have  moment  of  inertias  /„  =  15,  I,=  18,  and  /,=  3,  with 
K,  =  Kf  —  2.0,  K,  =  1.0,  and  r,  =  xy  =  r,  =  1.0.  Figure  33  on  page  45  shows  <f>  (position 
error),  <f>  (velocity),  and  the  restoring  torque  (T^J  generated  by  the  reaction  wheel  with 
respect  to  time.  Figures  34  and  35  show  8,  6  ,  T^,  and  ^  ,  T^,  respectively.  With 
the  relatively  small  moments  of  inertias,  stabilization  of  all  axes  within  the  ±  1°  criteria 
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occurs  within  40  seconds.  Increasing  the  moments  of  inertias  to  /,  =  95,  lf  =  98  ,  and 
/,  =  3  while  holding  A'„  Ky  at  2.0,  A',  at  1.0,  and  t„  t,  ,  and  t2  constant  at  1.0,  yields 
Figures  36-3S.  As  expected,  the  </>  and  0  axes  take  longer  to  stabilized  (approximately 
160  seconds)  than  the  above  case  with  the  smaller  moments  of  inertia.  The  t axis  of 
Figure  38  stabilizes  much  faster  than  either  the  <£  or  0  axes.  This  is  due  to  the  relatively 
small  moment  of  inertia  (/,)  compared  to  the  restoring  torque  generated  by  the  reaction 
wheel. 


PROGIIAM/UHI T  HAH! 

. 

SP  Reaction  Vhetl 

unit  Ttri  no.  (Mr) 

; 

369000  j 

MOMENTUM  (FT -LB -SEC) 

t 

o.A  t  1250  RrM 

MIGHT  (LOS.) 

f 

9.9 

sun 

: 

6.0"  DIA.  MAI  1  3.9"  RICH  (excluding  Mounting 
feet  projection) 

STALL  Tonqus  (OI-1N) 

1 

7,7 

TEA A  TOMER  (MATTS) 

! 

6.5 

MAX  STEED  (RTN) 

I 

1350 

MOTOR  EXCITATION 

I6V,  0-feak,  Squartvaie,  9oo  lit,  2  Thaae 

STEADT  STATE  TOMER(  MATTS ) ! 

2.5  1  1350  RTM 

TACHOMETER  TTTE 

T 

r.N.  Rotor,  2  Thaae  Mound  Stator,  Slnuaolda) 

Output,  I  CicWRo'/Ttiaat ,  Amplitude  warlea 
with  apeed. 

QUAHTITT  BUILT 

: 

92 

COMMENTS 

! 

AC  Squirrel  Cage  Induction  Motor 

FLIGHT  III  STORT 

More  than  20  oT  theae  unlta  hare  flown. 

Un\ta  with  operating  llae  treater  than  S 
jeara  are  preaentljr  functioning  ln-orblt. 

Figure  32.  BendLx  Aerospace  Reaction  Wheel  Technical  Data. 
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Figure  34.  Gravity  gradient  effect  with  3  RYV  (Y  axis)  Ix=l5,  Iy=  18,  Iz=3 
K,=  2.0  ,  Tj  =  1.0 
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Figure  38.  Gravity  gradient  effect  with  3  RW  (Z  axis)  Ix=95,  Iy  =  98,  Iz=3 
K,=  1.0  ,  t,=  1.0 
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Figures  39-41  have  moments  of  inertia  of  /,  =  15,  7,  =  18,  and  I,  =  3  with 
K,  =  Ky  =  .2,  K2  =  .1,  and  r,  =  r,  =  t,  =  1.0.  The  reduction  of  A' ,  A,,  and  A*  reduces  the 
reaction  wheel  torque  generated. 

Figure  39  shows  ^(position),  <£( velocity),  and  reaction  wheel  torque  (Tmt) 
with  respect  to  time.  Figures  40  and  41  show  6,  6,  T„,e  ,  and  \j/,  ,  respectively. 

Stabilization,  within  the  ±1°  criteria,  occurs  for  all  axes  in  under  200  seconds.  Com¬ 
parison  of  the  magnitudes  of  the  torque  values  of  Figures  33-35  and  Figures  39-41  for 
the  same  moment  of  inertia  values  shows  a  decrease  of  roughly  a  factor  of  five. 

Increasing  the  moments  of  inertia  to  I ,  =  95,  Iy  =  98  and  I2  =  3,  while  hold¬ 
ing  Kx=  Ky  =  .2,  K,  —  .1,  and  r,  =  ry  =  t,  =  1.0,  produces  Figures  42-44.  This  last  group 
of  reaction  wheel  Figures  shows  that  «650  seconds  are  required  to  dampen  the  system 
to  the  ±  1°  stabilization  criteria.  The  tradeoffs  between  the  different  reaction  wheel 
combinations  are  essentially  time  response  against  torque/power.  The  faster  the  desired 
response,  the  greater  the  magnitude  of  the  restoring  torque,  and  ultimately,  the  more 
power  that  will  be  consumed  in  achieving  this  response. 
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Figure  40.  Gravity  gradient  effect  with  3  RW  (Y  axis)  Ix=  15,  Iy=  18,  lz=  3 
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Gravity  gradienT^T^  — - ___ 
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d.  Effects  of  Magnetic  torquers 

For  the  magnetic  torquing  simulations,  the  initial  conditions  are  0,  8,  and 
<//  position  errors  =  5°  ,  and  <£,  6,  and  i//  =  0.  The  initial  restoring  torque  is  provided 
by  the  gravity  gradient  effect  combined  with  the  rigid  body  dynamics.  For  the  following 
simulations,  one  of  the  two  conditions  must  exist  prior  to  energizing  the  magnetic  coils 
to  generate  magnetic  torque: 

1)  position  greater  than  .5°  and  velocity  greater  than  0,  or 

2)  position  less  than  —.5°  and  velocity  less  than  0.  For  the  following  figures, 
the  magnetic  moment  varies  from  9.0427  to  .5  Wb-m.  At  the  maximum  value  of 
9.0427W-m,  N  =  400  turns,  I  =  .5  amp,  and  A  =  nr 1  where  r  is  the  radius  at  .  1 2m.  Subse¬ 
quent  values  are  obtained  by  reducing  the  current  drawn. 

The  magnetic  field  in  component  form  for  an  equatorial  orbit  can  be  ap¬ 
proximated  by  eqs(  110-11 2).  Since  the  earth's  magnetic  field  is  highly  dependent  on  the 
satellite's  position.  ORION"  must  have  the  capability  to  determine  its  position  in  order 
to  effectively  utilize  magnetic  torquing.  For  simulation  purposes,  the  satellite  is  assumed 
to  be  directly  over  Greenwich  at  time  t  =  0.  The  magnetic  coils  are  turned  on  and  off  at 
the  specified  strength  as  needed.  Varying  the  strengths  of  the  magnetic  moments  gener¬ 
ated  would  enhance  performance,  but  is  beyond  the  scope  of  this  thesis.  The  assumption 
that  magnetic  torque  may  be  generated  in  any  direction  due  to  vector  addition  is  not 
always  valid:  magnetic  torque  may  not  be  generated  in  the  direction  of  the  magnetic 
field. 

Figures  45-47  utilized  magnetic  moments  of  mx  —  my  =  mz  —  9.0427  with 
moments  of  inertias  IM  =  15,  Iy  =  18,  and  I ,  =  3.  For  Figure  45  on  page  60,  the  gravity 
gradient  restoring  torque  takes  approximately  800  seconds  before  meeting  a  condition 
necessary  for  energizing  the  magnetic  coils.  Due  to  the  small  moments  of  inertia,  the  <t> 
restoring  torque  is  needed  for  only  a  short  time.  Figures  46  and  47  show  the  6  and  <// 
axes  essentially  stabilized  within  the  ±1°  criteria  after  the  the  initial  magnetic  torque 
pulse.  For  the  8  axis,  this  pulse  occurs  at  the  1200  seconds  and  lasts  for  approximately 
300  seconds.  The  8  axis  has  a  lower  frequency  response  and,  therefore,  takes  longer  to 
to  meet  the  condition  necessary  to  energize  the  magnetic  coils.  Due  to  its  relatively  small 
moment  of  inertia,  the  axis  chatters  slightly  before  stabilizing.  A  smaller  magnitude 
of  restoring  torque  eliminates  this  problem. 

Increasing  the  moments  of  inertia  to  I ,  =  95,  /,  =  98,  and  4  =  3  while  leav¬ 
ing  the  magnetic  moments  mx,  my,  and  m,  constant  at  9.0427wb-m,  produces  Figures 
48-50.  As  expected,  more  restoring  torque  is  needed  to  stabilized  the  <j>  and  8  axes  forc- 
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ing  the  magnetic  coils  to  remain  energized  for  a  longer  period  of  time.  The  time  needed 
for  all  axes  to  be  within  the  ±  1°  criteria  is  limited  by  the  9  axis.  Stabilization  occurs 
within  2S00  seconds  while  for  the  case  with  the  smaller  moments  of  inertias,  stabilization 
occurs  at  approximately  2100  seconds. 
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Figure  45.  Gravity  gradient  with  magnetic  torquing  (X  axis):  Magnetic 


moment  =  9.0427W-m  lx=  15,  ly  =  18,  Iz=  3. 
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Figure  46.  Gravity  gradient  with  magnetic  torquing  (Y  axis):  Magnetic 
moment  =  9.0427W-m  ix=  15,  ly=  18,  Iz=3. 
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Figure  47.  Gravity  gradient  with  magnetic  torquing  (Z  axis):  Magnetic 
moment  =  9.0427W-m  Ix=  15,  ly  =  18,  Iz=  3. 
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Figure  48.  Gravity  gradient  with  magnetic  torquing  (X  axis):  Magnetic 
moment  =  9.0427W-m  Ix=  95,  Iy  =  98,  Iz=  3. 
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Figure  49.  Gravity  gradient  with  magnetic  torquing  (Y  axis):  Magnetic 
moment  =  9.0427W-m  lx  =  95,  ly  =  98,  lz  =  3. 
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Figure  50. 


Gravity  gradient  with  magnetic  torquing  (Z  axis):  Magnetic 
moment  =  9.0427W-m  Ix=95,  I y  =  98,  Iz=3. 
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Reducing  the  strength  of  the  magnetic  moments  reduces  the  torque  that  is 
generated  and.  thereby,  slows  the  response  time  to  stabilization.  For  Figures  51-56,  the 
magnetic  moments  m ,  =  my  =  ml  =  1.0.  Figures  51-53  have  moment  of  inertias  l,  =  15, 
4  =  18  ,  and  4  =  3.  Stabilization  for  the  6  axis  occurs  within  approximately  3600  sec¬ 
onds.  Increasing  the  moments  of  inertias  to  4  =  95,  Jy  =  98,  and  4  =  3  produces  Figures 
54-56.  Stabilization  within  the  +  1°  criteria  for  the  d  axis  occurs  at  approximately  11000 
seconds.  The  length  of  time  that  the  magnetic  torqers  are  energized  increases  as  the 
strength  of  the  moment  decreases.  Figures  57-62  show  the  effects  of  magnetic  torquing 
using  a  maximum  magnetic  moment  of  .5Wb-m  for  each  coil.  Figures  57-59  show  that 
for  low  moments  of  inertia  (4  =  15,  Iy  =  18,  and  4  =  3),  stabilization  occurs  within  the 
±  1°  within  48C0  seconds.  Increasing  the  moments  of  inertia  to  4  =  95,  4  =  98,  and 
4  =  3  holding  the  maximum  magnetic  moment  =.5Wb-m,  produces  Figures  60-62.  In 
Figure  61,  the  9  axis  may  take  a  prohibitively  long  time. 
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Figure  54.  Gravity  gradient  with  magnetic  torquing  (X  axis):  Magnetic 
moment  =  l.OW-m  Ix  =  95,  ly  =  98,  Iz=  3. 
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Figure  57. 


Gravity  gradient  with  magnetic  torquing  (X  axis):  Magnetic 
moment=  .50W-m  lx  =  15,  Iy  =  18,  Iz=  3. 
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Figure  59.  Gravity  gradient  with  magnetic  torquing  (Z  axis):  Magnetic 
moment  =  .5U\V-m  Ix=  15,  ly  =  18,  Iz=  3. 
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Figure  60.  Gravity  gradient  with  magnetic  torquing  (X  AXIS):  Magnetic 
moment  =  .50W-m  lx  =  95,  Iy  =  98,  Iz  =  3. 


Figure  61.  Gravity  gradient  with  magnetic  torquing  (Y  AXIS):  Magnetic 

moment  =  .50W-m  ix  =  95,  Iy  =  98,  Iz  =  3. 


77 


3.  Gravity  Gradient  capture  (from  tumbling  mode) 

Increasing  the  inertia  wiil  reduce  the  tumbling  rate  by  conforming  to  the  con¬ 
servation  of  angular  momentum,  /0u)0  =  /,&>,.  As  the  inertia  increases,  a  corresponding 
decrease  in  angular  rate  must  be  realized.  This  is  useful  in  obtaining  the  initial  gravity- 
gradient  capture  from  a  tumbling  mode.  Figure  63  on  page  80  shows  that  for  a  slow 
tumbling  satellite,  the  change  in  the  moment  of  inertia  by  extending  the  boom  may  be 
sufficient  to  allow  gravity  gradient  capture.  The  satellite  is  initially  tumbling  at  a  rate 
of  5  rotations  per  orbit  about  the  pitch  axis.  The  extension  of  the  boom  increases  the 
moment  of  inertia  about  the  y  axis  by  a  factor  of  6;  therefore,  by  the  conservation  of 
angular  momentum,  the  angular  rate  must  decrease  by  a  factor  of  6.  This  reduces  the 
rotations  per  orbit  to  approximately  5  6  per  orbit  well  below  the  one  rotation  per  orbit 
rate  necessary  to  acheive  gravity  gradient  capture.  Figure  63a  shows  the  initial  tumbling 
rate  of  5  rotations  per  orbit.  Figure  63b  shows  the  corresponding  decrease  in  the  tum¬ 
bling  rate  as  the  moment  of  inertia  is  increased.  Figure  63c  shows  the  increase  in  mo¬ 
ment  of  inertia  as  the  boom  is  extended  with  respect  to  time. 

When  an  initial  tumbling  rate  is  higher  than  the  rate  that  can  be  captured  solely 
through  the  extension  of  the  boom,  magnetic  torquing  may  aid  in  gravity  gradient  cap¬ 
ture.  By  creating  torque  in  the  direction  opposite  to  that  of  the  angular  velocity,  mag¬ 
netic  torquing  can  be  used  to  slow  the  rate  of  the  tumbling  vehicle. 
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Figure  63.  Gravity  gradient  capture  by  extension  of  boom 


V.  CHAPTER  5 


A.  CONCLUSIONS 

1.  Gravity  gradient  stabilization  (only) 

Gravity  gradient  stabilization  alone  is  not  a  viable  attitude  control  method  for 
ORION.  This  is  due  to  gravity  gradient  effect  providing  little  yaw  restoring  torque.  In 
addition,  gravity  gradient  stabilization  provides  little  dampening  effect  on  the  roll  and 
pitch  axes.  When  augmented  by  other  control  schemes,  gravity  gradient  stabilization 
provides  sufficient  passive  restoring  torque  to  warrant  its  use. 

2.  Gravity  gradient  with  3  reaction  wheels 

a.  Advantages 

Gravity  gradient  stabilization  augmented  by  three  reaction  wheels,  each 
aligned  along  the  body's  axes,  provides  excellent  attitude  control  results  both  in  time 
response  and  accuracy. 

b.  Disadvantages 

Power  consumption  may  be  prohibitively  high  for  this  method.  At  steady 
state,  each  reaction  may  use  up  to  6  watts  while  drawing  considerably  more  during  the 
acceleration  or  deceleration  phase.  If  a  steady  state  drain  of  18  watts  (6  for  each  re¬ 
action  wheel)  is  acceptable,  then  this  method  should  be  considered.  The  projected  end- 
of-life  power  budget  for  ORION  is  60  watts  [Ref.  1,  p.  132].  Using  this  figure  as  an 
estimate,  a  steady  state  draw  of  18  watts  leaves  42  watts  available  for  mission  support. 
From  Figure  3  on  page  3,  an  average  of  34  watts  is  required  for  mission  support.  The 
size  and  weight  of  three  reaction  wheels  may  be  prohibitively  high  as  well.  Another 
disadvantage  is  that  thrusters  must  be  used  to  remove  momentum  when  the  reaction 
wheel(s)  become  saturated, 

3.  Control  moment  gyro  (CMG) 

a.  Advantages 

Al^ough  not  discussed  in  detail,  control  moment  gyros  are  capable  of 
generating  large  restoring  torques. 

b.  Disadvantages 

The  weight  requirements  for  a  control  moment  gyro  are  prohibitively  high 
when  compared  to  reaction  wheels  when  relatively  small  torques  are  invoved.  As  a  small 
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multi-purpose  satellite,  ORION  does  not  need  the  magnitude  of  restoring  torque  neces¬ 
sary  to  warrant  the  use  of  control  moment  gyros. 

4.  Gravity  gradient  with  three  magnetic  torquers 

a.  Advantages 

Three  othogonally  oriented  magnetic  torquers  have  the  advantage  of  gen¬ 
erating  restoring  torque  in  almost  any  direction  based  on  vector  addition.  In  this  man¬ 
ner,  the  corrective  torque  may  be  produced  to  counter  any  disturbance  from  any 
direction  as  long  as  it  is  small  enough  not  to  tumble  the  vehicle.  Another  advantage  for 
magnetic  torquing  is  that  the  torquers  only  draw  power  while  providing  restoring  torque. 
This  amounts  to  a  considerable  power  savings  when  compared  to  the  reaction  wheel 
The  power  drawn  by  the  magnetic  torquers  during  operations  is  based  on  the  amperes 
used.  This  allows  for  flexibility  of  design  since  doubling  the  number  of  turns  and  halving 
the  current  drawn  still  provides  the  same  magnitude  of  the  magnetic  moment. 

b.  Disadvantages 

Magnetic  torquers  are  altitude  dependent.  The  magnitude  of  the  earth's 
magnetic  field  is  a  function  of  -jy-.  At  high  altitudes,  magnetic  torquing  is  not  effective. 

A* 

In  addition,  magnetic  torquing  is  not  always  available  in  all  directions.  Recalling  that 
magnetic  torque  is  Tmei  =  M  x  B,  magnetic  torque  may  not  be  generated  in  the  direction 
of  the  B  field  due  to  vector  cross  product  relationships.  Due  to  the  continuously  chang¬ 
ing  B  field  with  respect  to  the  satellite,  this  is  a  minor  disadvantage.  The  desired  restoring 
torque  can  soon  be  generated  as  soon  as  the  B  field  direction  changes. 

5.  Thrusters 

a.  Advantages 

At  high  altitudes,  thrusters  and  momentum  exchange  devices  are  the  only 
control  methods  available  for  accurate  stabilization.  In  addition,  thrusters  are  the  most 
effective  method  for  stabilizing  a  vehicle  with  a  high  tumbling  rate. 

b.  Disadvantages 

Thrusters  are  necessary'  for  satellites;  however,  any  fuel  used  for  attitude 
control  or  momentum  removal  is  less  fuel  available  for  orbit  transfer  and  maintainance. 
For  a  small,  lightweight  satellite  with  limited  fuel  storage  capacity,  fuel  resources  must 
be  economized. 

6.  Overview- 

In  the  final  analysis,  the  attitude  control  system  employed  will  be  based  on 
specific  mission  requirements.  Only  when  the  ORION  mission  is  known  can  proper 
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weights  be  given  to  each  of  the  attitude  control  system  design  factors,  most  noticeably, 
power  consumption,  pointing  accuracy,  altitude,  and  desired  lifetime. 
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APPENDIX  GRADIENT  STABILIZATION  PROGRAM 


TITLE  3  AXIS  STABILIZATION  WITH  GRAVITY  GRADIENT  EFFECT 


*  DEFINITION  OF  TERMS 

* 


*  ALT=ALTITUDE  OF  SATELLITE  (METERS) 

*  R1=RADIUS  OF  EARTH  (METERS) 

*  R=RADIUS  OF  EARTH  +  SATELLITE  ALT  (METERS) 

*  RO=DENSITY  OF  AIR  AT  THE  ORBITAL  ALTITUDE 

*  W0=ANGULAR  RATE  OF  SATELLITE  WITH  RESPECT  TO  EARTH 

*  D=DRAG  FORCE=l. 7E+14*(RO/R) 

*  PS I, THETA, PHI  ARE  ANGLES  ABOUT  THE  X,Y,Z,  AXES  (RADIANS) 

*  PSIC ,THETAC ,PHIC  ARE  ANGLES  ABOUT  THE  X.Y.Z  AXES  (DEGREES) 

*  PSID,TD,PHID  ARE  ANGULAR  VELOCITIES  ABOUT  THE  X,Y,Z  AXES  (RADS/SEC) 

*  (Bl) ,B(2) ,B(3)  ARE  THE  ANGULAR  ACCELERATION  ABOUT  THE 

*  X,Y,Z  AXES  (RADS/SEC**2) 

*  IX, IY, IZ  ARE  THE  MOMENTS  OF  INERTIA  ABOUT  THE  X,Y,Z  AXES  (KG-M**2) 

*  TX,TY,TZ  ARE  THE  SUMMATION  OF  THE  DISTURBANCE  TORQUES  AND 

*  RESTORING  TORQUES  EXCLUDING  GRAVITY  GRADIENT  EFFECTS 

*  TMX,TMY,TMZ  ARE  MAGNETIC  TORQUES  ABOUT  THE  X,Y,Z  AXES 

*  TRWX , TRWY , TRW Z  ARE  REACTION  WHEEL  TORQUES  ABOUT  THE  X, Y,Z  AXES 

*  TSX,TSY,TSZ  ARE  SOLAR  PRESSURE  TORQUES  ABOUT  THE  X,Y,Z  AXES 

*  TAX,TAY,TAZ  ARE  AERODYNAMIC  TORQUES  ABOUT  THE  X,Y,Z  AXES 

*  LX , LY , LZ  IS  THE  DISTANCE  BETWEEN  THE  CENTER  OF  PRESSURE  AND 

*  CENTER  OF  MASS  (M) 

*  MX,MY,MZ=SPACECRAFT’S  GENERATED  MAGNETIC  MOMENT  WITH  RESPECT  TO 

*  THE  X,Y,Z  AXES  (WB-M)  M=NIA  WHERE  M=MAGNETIC  MOMENT, 

*  N=#  OF  TURNS,  I=CURRENT  (AMPS),  AND  A=AREA  (M**2).  FOR 

*  THE  MAXIMUM  VALUE  OF  M=9.0427,  N=400 , 1=  5A,RADIUS=.  12M. 

*  MG-STRENGTH  OF  EARTH'S  MAGNETIC  DIPOLE=8.  0E15  WB-M  AT  R=0 

*  BX , BY , BZ=EARTH ' S  MAG  FIELD  WITH  RESPECT  TO  X,Y,Z(WB/M**2) 

*  LAMS=RIGHT  ASCENSION  FROM  GREENWICH  @  TIME=0,  IN  DEGREES( LONGITUDE) 


ARRAY  B( 3) 

CONST  ALT=1000E+03 ,RO=7.  0E-10 ,R1=6378. E+03,.  .  . 

IX=95.  ,  IY=98.  ,  IZ=03. 00,  LAMS0=0,... 

MX=9.  0478 ,  MY=9.  0427 ,  MZ=9.0427,  MG=8.  OE+15 , .  .  . 

KX=2.  ,  KY=2.  ,  KZ=1.  ,  TAU1=1.  ,  TAU2=1.  ,  TAU3=1 . 

TSX=0.  ,  TSY=0. ,  TSZ=0.  ,  TAX=0.  ,  TAY=0. ,  TAZ=0. 

* 

*  PHI, THETA, AND  PSI  ARE  INIT  5  DEGREES  =  .08725  RADIANS  OF  ERROR* 

********************************************************************** 


INCON  ITHETA=.  08725, IPHI=.  08725 ,IPSI=.  08725 , ILAMS=0 ,PHIC=5.  . 
THETAC=5.  , PSIC=5. 

* 


DERIVATIVE 

R=R1+ALT 

D=l. 7E+14*(R0/R) 
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WO=SQRT( 3.  987E+14/(R**3) ) 


NOSORT 

LAMS=LAMSO+WO*TIME*( 57.  3) 
10  IF  (LAMS.  GE.  360)  THEN 
LAMS=LAMS-360. 

GO  TO  10 
ENDIF 


******************************************************* 

*  THIS  SECTION  ADDS  THE  REACTION  WHEEL  EQUATIONS  * 

******************************************************* 

*  TRWX= - KX* ( TAU 1*PH I D  +  PHI) 

*  TRWY= - KY* ( TAU2*TD  +  THETA) 

*  TRWZ=-KZ*(TAU3*PSID  +  PS I) 

***********£ nd  REACTION  WHEEL************************** 
******************************************************* 

*  THIS  SECTION  ADDS  THE  MAGNETIC  TORQUE  EQUATIONS  * 

*  AND  SOLVES  AN  EXAMPLE  FOR  THE  EARTH'S  MAGNETIC  * 

*  FIELD  (EQUATORIAL  ORBIT  ASSUMED)  * 

*iV********“VVfV?****Vr,jV*'»V,>VVrVrVf’>VVc^r,>V'jVVf,>VVfiV>V'»V‘>V,3V*7V^r,A'*3Sp,>VVr^V^r>V'iV’jVil'rVt,VrTV?V 

BX=(MG/(R**3) )*SIN( .  2042)*SIN( ( ( LAMS*PI)/ 180) -5.  0789) 
BY=-(MG/(R**3))*COS(.  2042) 

BZ=-(MG/(R**3))*SIN(. 2042)*COS(((LAMS*PI)/180)-5. 0789) 


TMAGX=( MZ*B Y ) - ( MY*B  Z ) 

TMAGY=( MX*BZ) - ( MZ*BX) 
TMAGZ=(MY*BX)-(MX*BY) 

****************************** 

*  CONDITIONS  FOR  TORQUE  GENERATION  MET? 

***** ************************* 

IF( PHIC.  GT.  .  5)  THEN 
IF( PHID. GT. 0)  THEN 
TMX=-( ABS(TMAGX) ) 

ELSE 

TMX=0 

ENDIF 

ENDIF 

IF(PHIC.  LT.  -.5)  THEN 
IF( PHID. LT.  0)  THEN 
TMX=AB  S ( TMAGX ) 

ELSE 

TMX=0 

ENDIF 

ENDIF 

IF(THETAC.  GT.  .  5)  THEN 
IF(TD. GT. 0)  THEN 
TMY=-( ABS(TMAGY) ) 

ELSE 

TMY=0 

ENDIF 

ENDIF 

IF(THETAC. LT. -.5)  THEN 
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IF(TD.  LT.  0)  THEN 
TMY=AB  S ( TMAGY ) 

ELSE 
TMY=0 
END  IF 
END  IF 

IF(PSIC.GT.  .  5)  THEN 
IF(PSID.  GT.  0)  THEN 
TMZ=-(ABS(TMAGZ) ) 

ELSE 
TMZ=0 
END  IF 
END  IF 

IF( PSIC. LT.  5)  THEN 
IF(PSID.  LT.  0)  THEN 
TMZ=ABS ( TMAGZ ) 

ELSE 
TMZ=0 
END  IF 
ENDIF 

*********************END  MAGNETIC  TORQING  SECTION******************** 


TX=TMX+TRWX+TSX+TAX 

TY=TMY+TRWY+TSY+TAY 

TZ=TMZ+TRWZ+TSZ+TAZ 


*******F OR  GRAVITY  GRADIENT  EFFECT  ONLY  SET  TX=TY=TZ=0************** 


B(1)=-4*(W0**2)*((IY-IX)/IZ)*PHI+W0*((IX+IZ-IY)/IX)*PSID+(TX/IX) 
B( 2)=-3*( W0**2)*( ( IX-IZ) /IY)*THETA+TY/IY 

B(3)=-(W0**2)*((IY-IX)/IZ)*PSI-W0*((IX+IZ-IY)/IZ)*PHID+(TZ/IZ) 

* 

* 

PHID=INTGRL( 0.  ,B(1)) 

TD=INTGRL( 0. ,B(2)) 

PSID=INTGRL( 0. ,B(3)) 


CPHI=INTGRL( 0.  ,PHID) 

CTHETA=INTGRL( 0.  ,TD) 

CPSI=INTGRL( 0.  ,PSID) 

* 

PHI=IPHI+CPHI 

THETA=ITHETA+CTHETA 

PSI=IPSI+CPSI 

*PH I, THETA, PS I  ARE  =  TO  INITIAL  CONDITIONS  AS  ANNOTATED  WITH  AN  "I” 
*PRECEDING  VARIABLE+CALCULATED  VALUES  AS  ANNOTATED  WITH  THE  "C"  IN  FRONT 


PHIC=PHI*57.  3 
THETAC=THETA*5  7 .  3 
PSIC=PSI*57. 3 

* 

***************************+** ************************************ 


*  USE  FOLLOWING  IF  GRAVITY  GRADIENT  EFFECT  ONLY 
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*****************************************  ************************* 
*ONTRL  FINTIM=20000 ,DELT=200. 

*RINT  200.  ,  PHIC,  THETAC,  PSIC,  PHIDC,  TDC,  PSIDC 
*AVE  200.,  PHIC,  THETAC,  PSIC,  PHIDC,  TDC,  PSIDC 
*RAPH  (G1 ,DE=TEK618)  TIME(UN=SEC) ,PHIC(UN=DEGREES) 

*ABEL  (G1,DE=TEK618)  PHI(DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G2 ,DE=TEK618)  TIME ( UN=SEC ) , THETAC ( UN=DEGREES ) 

*ABEL  (G2,DE=TEK618)  THETA(DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G3,DE=TEK618)  TIME(UN=SEC) ,PSIC(UN=DEGREES) 

*ABEL  (G3 ,DE=TEK618)  PSI(DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G4,DE=TEK618)  TIME(UN=SEC) ,PHIDC(UN=' DEG/SEC ' ) 

*ABEL  (G4,DE=TEK618)  PHID  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G5,DE=TEK618)  TIME(UN=SEC) ,TDC(UN=' DEG/SEC' ) 

*ABEL  (G5,DE=TEK618)  TDC  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3.  00 
*RAPH  (G6,DE=TEK618)  TIME(UN=SEC) ,PSIDC(UN=' DEG/SEC' ) 

*ABEL  (G6,DE=TEK618)  PSIDC  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3.  00 
*************END  GRAVITY  GRADIENT  ONLY  GRAPHS  ******************* 

Vc 

*  USE  THE  FOLLOWING  IF  REACTION  WHEELS  ARE  USED  * 

*********************************************************** 

*ONTRL  FINTIM=400,DELT=5. 

*RINT  100.  ,  PHIC,  THETAC,  PSIC,  TX,  TY,  TZ 

*AVE  5.  ,  PHIC,  THETAC,  PSIC,  TX,  TY,  TZ,  PHID,  TD,  PS ID 

*RAPH  (G1,DE=TEK618)  TIME(UN=SEC) , PHIC(UN=DEGREES) 

*ABEL  (G1,DE=TEK618)  PHI (DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G2,DE=TEK618)  TIME(UN=SEC) ,THETAC(UN=DEGREES) 

*ABEL  (G2,DE=TEK618)  THETA(DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G3,DE=TEK618)  TIME(UN=SEC) ,PSIC(UN=DEGREES) 

*ABEL  (G3 ,DE=TEK618)  PSI(DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G4,DE=TEK618)  TIME(UN=SEC) ,TX(UN=NM) 

*ABEL  (G4,DE=TEK618)  TXRW  VS  TIME  (GG<ai000K)  IX=95  IY=98  IZ=3.  00 
*RAPH  (G5,DE=TEK618)  TIME(UN=SEC) ,TY(UN=NM) 

*ABEL  (G5,DE=TEK618)  TYRW  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G6,DE=TEK618)  TIME(UN=SEC) ,TZ(UN=NM) 

*ABEL  (G6,DE=TEK618)  TZRW  VS  TIME  (GG<aiOOOK)  IX=95  IY=98  IZ=3.  00 
*RAPH  (G7,DE=TEK618)  TIME(UN=SEC) ,PHID(UN=' RAD/SEC ' ) 

*ABEL  (G7 ,DE=TEK618)  PHID  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G8,DE=TEK618)  TIME(UN=SEC) ,TD(UN=’RAD/SEC’ ) 

*ABEL  (G8,DE=TEK618)  TD  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G9 ,DE=TEK618)  TIME(UN=SEC) ,PSID(UN=' RAD/SEC' ) 

*ABEL  (G9,DE=TEK618)  PSID  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*********************£Nj)  REACTION  WHEEL  GRAPHS********************** 

*  * 
******************************************************************** 

*  USE  FOLLOWING  GRAPHS  IF  MAGNETIC  TORQUING  IS  UTILIZED  * 

******************************************************************** 

*ONTRL  FINTIM=4000.  ,DELT=5. 

*RINT  100.,  PHIC,  THETAC,  PSIC,  WO ,B( 1) ,B( 2) ,B( 3) , LAMS 

*AVE  5.,  PHIC,  THETAC,  PSIC,  LAMS,  TX, TY,TZ,BX, BY, BZ, PHID, TD, PSID 

*RAPH  ( G 1 ,DE=TEK618)  TIME(UN=SEC) ,PHIC(UN=DEGREES) 

*ABEL  (G1,DE=TEK618)  PHI(DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G2,DE=TEK618)  TIME(UN=SEC) ,THETAC(UN=DEGREES) 

*ABEL  (G2,DE=TEK618)  THETA(DEG)  VS  TIME  (GG(ai000K)  IX=95  IY=98  IZ=3.  00 
*RAPH  (G3,DE=TEK618)  TIME(UN=SEC) , PSIC(UN=DEGREES) 

*ABEL  (G3,DE=TEK618)  PSI(DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
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*RAPH  (G4,DE=TEK618)  TIME(UN=SEC) ,LAMS(UN=DEGREES) 

*ABEL  (G4,DE=TEK618)  LAMS  (DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3.  00 
*RAPH  (G5 ,DE=TEK618)  TIME(UN=SEC) ,TX(UN=N-M) 

*ABEL  (G5 ,DE=TEK618)  TX  VS  TIME  (GG@1000K) 

*FAPH  (G6,DE=TEK618)  TIME(UN=SEC) ,TY(UN=N-M) 

*ABEL  (G6 ,DE=TEK618)  TY  VS  TIME  (GG@1000K) 

*RAPH  (G7 ,DE=TEK618)  TIME(UN=SEC) ,TZ(UN=N-M) 

*ABEL  (G7,DE=TEK618)  TZ  VS  TIME  (GG@1000K) 

*RAPH  (G8,DE=TEK618)  TIME(UN=SEC) ,BX(UN=W-M) 

*ABEL  (G8,DE=TEK618)  BX  VS  TIME  (GG@1000K) 

*RAPH  (G9 ,DE=TEK618)  TIM£(UN=SEC) ,BY(UN=W-M) 

*ABEL  (G9,DE=TEK618)  BY  VS  TIME  (GG@1000K) 

*RAPH  (G10,DE=TEK618)  TIME(UN=SEC) ,BZ(UN=W-M) 

*ABEL  (G10,DE=TEK618)  BZ  VS  TIME  (GG@1000K) 

*RAPH  (G11,DE=TEK618)  TIME(UN=SEC) ,PHID(UN=' RAD/SEC' ) 

*ABEL  (Gil ,DE=TEK618)  PHID(RADSEC)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3.  00 
*RAPH  (G12 ,DE=TEK618)  TIME(UN=SEC) ,TD(UN=' RAD/SEC1 ) 

*ABEL  (G12 ,DE=TEK618)  TD(RADSEC)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 
*RAPH  (G13,DE=TEK618)  TIM£(UN=SEC) ,PSID(UN=' RAD/SEC1 ) 

*ABEL  (G13 ,DE=TEK618)  PSID(RADSEC)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 

END 

STOP 
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